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Abstract 

In tliis paper we study tlie interplay between tlie recently proposed F-tlieory 
GUTs and cosmology. Despite tlie fact that the parameter range for F-theory GUT 
models is very narrow, we find that F-theory GUTs beautifully satisfy most cos- 
mological constraints without any further restrictions. The viability of the scenario 
hinges on the interplay between various components of the axion supermultiplet, 
which in F-theory GUTs is also responsible for breaking supersymmetry. In these 
models, the gravitino is the LSP and develops a mass by eating the axino mode. The 
radial component of the axion supermultiplet known as the saxion typically begins 
to oscillate in the early Universe, eventually coming to dominate the energy density. 
Its decay reheats the Universe to a temperature of ~ 1 GeV, igniting BBN and 
diluting all thermal relics such as the gravitino by a factor of ~ 10"^ — 10^^ such 
that gravitinos contribute a sizable component of the dark matter. In certain cases, 
non-thermally produced relics such as the axion, or gravitinos generated from the de- 
cay of the saxion can also contribute to the abundance of dark matter. Remarkably 
enough, this cosmological scenario turns out to be independent of the initial reheat- 
ing temperature of the Universe. This is due to the fact that the initial oscillation 
temperature of the saxion coincides with the freeze out temperature for gravitinos in 
F-theory GUTs. We also find that saxion dilution is compatible with generating the 
desired baryon asymmetry from standard leptogenesis. Finally, the gravitino mass 
range in F-theory GUTs is 10 — 100 MeV, which interestingly coincides with the 
window of values required for the decay of the NLSP to solve the problem of ^Li 
over-production. 
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1 Introduction 



Two prominent triumphs of modern theoretical physics are the Standard Models of 
particle physics and cosmology. Moreover, the interplay between particle physics 
and astrophysics/early Universe cosmology has already proven to be a fruitful arena 
of investigation for both fields. On the astrophysics side, this interplay has led to 
the very successful predictions of Big Bang Nucleosynthesis (BBN), which accounts 
for the abundance of light elements. In addition, ideas from particle physics have 
provided several plausible mechanisms such as baryogenesis or leptogenesis which 
could generate the observed baryon asymmetry. Finally, many extensions of the 
Standard Model include dark matter candidates. 

On the particle physics side, constraints from astrophysics have led to novel, and 
sometimes quite stringent conditions on possible extensions of the Standard Model. 
These constraints can translate into important bounds on parameters of a candidate 
model which may be inaccessible from other avenues of investigation, and which can 
have repercussions beyond their immediate astrophysics applications. For example, 
compatibility with the successful predictions of BBN imposes important restrictions 
such as the requirement - spectacularly confirmed by LEP - that essentially there 
are at most three generations of light neutrinos! This is in amazing accord with 
the Standard Model of particle physics, and indicates a deep link between these 
two fields. Other cosmological considerations such as over-production of gravitinos 
in supersymmetric models, or a deficit in the observed baryon asymmetry provide 
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additional constraints. Satisfying all of these constraints is often a non-trivial task 
for a given model, but can also point the way to novel mechanisms which may not 
be available in the standard cosmology. 

At a more refined level, the interrelations between particle physics and cosmology 
roughly bifurcate into issues where gravity itself plays a central role, and questions 
where gravity plays only a supporting role in addressing more detailed features of a 
given particle physics model. For example, issues connected to the cosmological con- 
stant, or the homogeneity of the early Universe fall in the former category, whereas 
particle physics oriented issues such as the origin of dark matter or the overall baryon 
asymmetry fit most naturally in the latter category. Due to the vast array of obser- 
vational data from probes wholly separate from cosmology, issues more closely tied 
to particle physics appear to at present be more tractable. In this regard, it is there- 
fore quite natural, as is often done in the particle physics literature, to exclusively 
focus on the "low energy" aspects associated with the cosmology of a given particle 
physics model, parameterizing our ignorance of what occurs in the very early Uni- 
verse in terms of an "initial reheating temperature" for the Universe, T'^fj. For the 
purposes of particle physics considerations, the cosmology of the Universe effectively 
begins at this temperature. In many cases, this reheating temperature ends up being 
smaller than the Planck, or GUT scale, and models with as low as 10 — 100 TeV 
have also been discussed. 

It is interesting that this division in cosmology between gravity and particle 
physics issues parallels recent developments in the string theory literature. It is clear 
that a vast landscape of consistent string theory vacua exists. While this makes 
the string theory paradigm a very rich and flexible physical model, at present it 
also lacks predictive power because it does not lead to any particularly distinguished 
vacua! As a principle which can be used to limit the search for semi-realistic vacua, 
in [1-5] (see also [6-17]), it was shown that in the context of F-theory based mod- 
els, demanding the existence of a limit where the Planck scale (and thus associated 
gravitational questions) decouples, in tandem with some qualitative particle physics 
considerations such as the existence of a GUT, leads to a remarkably limited, and 
predictive framework. In fact, without adding any additional ingredients, tradition- 
ally vexing particle physics issues related to flavor hierarchies, the doublet-triplet 
splitting problem, the /i/i? /i-problem, and undesirable GUT mass relations all find 
natural solutions in F-theory GUT models. In addition, natural estimates for the 



3 



overall magnitude of the Yukawa couplings, axion decay constant, parameter and 
MSSM soft mass terms all fall in an acceptable phenomenological range. 

The aim of the present paper is to examine the cosmology of the F-theory GUT 
scenario as a model of particle physics. We find that these models naturally satisfy 
the typically stringent constraints derived from compatibility with BBN as well as 
bounds on the overall relic abundance of dark matter candidates, such as the grav- 
itino. Perhaps the most striking feature of this analysis is the absence of any major 
problem, namely, that consistency with cosmology imposes almost no constraint at 
all on the initial reheating temperature T^^. 

That this is the case is a highly non-trivial consequence of the parameter range 
found for F-theory GUTs. Indeed, one potentially significant source of tension could 
in principle have originated from the fact that as a model of gauge mediated super- 
symmetry breaking with a relatively high mass for the gravitino [3]: 

m3/2 = -^T^ ~ 10 - 100 MeV, (1) 

the relic abundance of gravitinos, which is the LSP, can potentially overdose the 
Universe. In the above, as throughout this paper, Mpi denotes the reduced Planck 
mass MpL = 2.4 x 10^^ GeV, and F denotes the component of the GUT group singlet 
chiral superfield X responsible for supersymmetry breaking: 

{X) = x + e^F (2) 

where as shown in [3], simultaneously solving the /i problem and generating viable 
soft mass terms in a minimal gauge mediation scenario requires: 

F ~ 10^^ GeV^ (3) 
X ~ 10^2 GeV. (4) 

We note that the Goldstino mode corresponding to the fermionic component of X is 
eaten by the gravitino. 

In many gauge mediation models, the scale of supersymmetry breaking is signif- 
icantly lower. From the perspective of cosmology, the relic abundance of gravitinos 
would at first appear to provide strong motivation for lowering the value of ^3/2. 
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Indeed, in the most straightforward approximation, it is quite natural to take as 
high as the GUT scale. In this case, a well known estimate for the relic abundance 
of gravitinos in many supersymmetric models requires: 



y»3/2 

2 keV 



< 0.1 



(5) 



to avoid an overabundance of gravitinos in the present Universe. In particular, this 
would appear to suggest an upper bound for F of order: 



which is significantly lower than equation ([3])! In the gravitino cosmology literature, 
it is common to take < Mqut 

~ 3 X 10^^ GeV to truncate the production of 
thermally produced gravitinos, thus avoiding precisely these types of issues. At a 
conceptual level, though, it is somewhat distressing that particle physics consider- 
ations demand such a stringent upper bound on T^^. Indeed, insofar as the value 
of T^jj is dictated by gravitational issues which are a priori wholly separate from 
details of a particular particle physics models, this type of tuning of parameters is 
quite puzzling. 

In F-theory GUTs, the resolution of this apparent dilemma again resides in the 
chiral superfield X. The essential point is that X plays a dual role in F-theory 
GUTs because its bosonic component breaks the anomalous global U{1) Peccei- 
Quinn symmetry of the low energy theory. As such, the associated Goldstone mode 
is the QCD axion, with decay constant: 



solving the strong CP problem. In the context of supersymmetric theories, however, 
the axion corresponds to one of two real degrees of freedom associated with the 
bosonic component of the corresponding axion supermultiplet. The other degree of 
freedom, known as the saxion is exactly massless in the limit where supersymmetry 
is restored, and in the present context has a mass of order 100 GeV. An exciting 
feature of the F-theory GUT is that the mass of the saxion is controlled by UV 
sensitive details of the compactification, such as the mass of the anomalous U{1)pq 
gauge boson. As such, cosmological constraints for the saxion provide a window into 



F < 7 X 10^^ GeV^ 



(6) 



fa ~ V2x ~ 10^2 GeV, 



(7) 
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the high scale dynamics of the model. 

Because the potential of the saxion is nearly flat, it is easily displaced from its 
minimum, and will generically begin to oscillate as the early Universe cools from the 
initial reheating temperature T^^j until the saxion decays. In a generic situation, the 
initial amplitude for sq is sufficiently large that its vacuum energy density comes to 
dominate the energy density of the Universe. The value of sq is on the order of the 
characteristic Kaluza-Klein scale for the X field: 

So ~ Mx ~ lO^^-^ GeV. (8) 

The decay of the saxion releases a large amount of entropy into the Universe, effec- 
tively diluting the overall relic abundance of all particle species. As a consequence, 
we find that rather neatly, one component of the axion supermultiplet, the saxion, 
counteracts the potentially dangerous features of the fermionic gravitino component 
(which includes the axino as the longitudinal degree of freedom)! We note that the 
decay of the saxion or some other cosmological modulus as a means to dilute the relic 
abundance of a particle species to acceptable levels has certainly been discussed in 
the literature before, for example in [18]. The primary novelty here is that without 
any additional assumptions, F-theory GUT models automatically resolve the most 
problematic features of gravitino cosmology. The requisite era of saxion dominance 
occurs provided the initial reheating temperature satisfies: 

> 10^ GeV. (9) 

In this case, the dilution of the saxion is such that the gravitino could naturally make 
up a prominent component of the dark matter density. Depending on whether the 
axion begins to oscillate before or after an era of saxion domination, the axion can 
also provide a component of the dark matter. We are currently investigating whether 
F-theory GUTs provide additional dark matter candidates [19]. Finally, we also find 
that lower values of the reheating temperature are also possible, and are compatible 
with a regime where the saxion has a smaller initial amplitude. In this case, the 
gravitino and axion can potentially both contribute to the dark matter density. 

But while the evolution of the saxion neatly solves the gravitino "problem" , it can 
in principle introduce additional constraints. For example, the decay of the saxion, 
with its significant entropy release might disrupt the start of BBN. Most conserva- 



6 



tively, this requires that the reheating temperature for the saxion remain above the 
starting temperature for BBN. Quite remarkably, this requirement is again satisfied 
quite naturally in F-theory GUT models. We find that the reheating temperature of 
the saxion, T^^j satisfies: 

T'^H ~ 1 GeV > Tbbn ~ 2 MeV. (10) 

Compatibility with BBN also imposes strong constraints on the decay products of 
the saxion. Using the well known bound (which we shall review) on the branching 
ratio of the saxion to relativistic species such as the axion, we find that we must 
either posit the existence of an additional decay product not found in the MSSM, or 
that the saxion must have sufficient mass to decay dominantly to two Higgs fields so 
that: 

TUsax > 2mho ~ 230 GeV. (11) 

A sufficient amount of baryon asymmetry must be generated at high tempera- 
tures in order for BBN to produce the observed abundances of light elements. In 
the standard solution to the gravitino problem, lowering the initial reheating tem- 
perature T]j^ has the deleterious consequence of removing some of the more efficient 
mechanisms for generating such an asymmetry, such as GUT scale baryogenesis, or 
standard leptogenesis. Rather, in this context it is common to invoke a mechanism 
where a large baryon asymmetry can be generated by the coherent oscillation of a 
field which carries either non-trivial baryon, or lepton number, such as in the Affleck- 
Dine scenario. In the present context, however, the dilution of the saxion completely 
relaxes any upper bound on T^^h- Thus, the only condition left to check is whether 
any of the available mechanisms are capable of generating a sufficiently large initial 
baryon asymmetry which can survive the effects of dilution. Performing this analysis, 
we find that standard leptogenesis indeed produces an appropriate baryon asymme- 
try. In other words, the range of parameters suggested by F-theory GUTs naturally 
fall within the small range of dilution factors compatible with diluting the graviti- 
nos to an amount where they do not overdose the Universe, but which nevertheless 
preserve a sufficient baryon asymmetry]^ 

The organization of this paper is as follows. In section [2] we review the main 



^See [20] for other recent work on moduli stabilization, leptogenesis, and dark matter in the 
context of gauge and gravity mediated scenarios. 
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interrelations between cosmology and particle physics which will occupy a central 
role in the present paper. In section [3] we study the axion supermultiplet and its 
interactions in the context of F-theory GUTs. Section H] forms the main body of our 
paper in which we present the particular cosmological scenario suggested by F-theory 
GUTs, paying special attention to the role of the saxion. In this same section, we 
analyze the gravitino and axion relic abundance, study constraints from BBN, and 
determine the overall baryon asymmetry generated by standard leptogenesis. Finally, 
in section [5] we briefly discuss future avenues of investigation in the cosmology of F- 
theory GUTs. 



2 Cosmology and Particle Physics 

One of the primary aims of the present paper is to determine how cosmological 
considerations constrain F-theory GUTs. In this section we review the main issues 
which must be addressed in a viable cosmological scenario. This section is entirely 
review, and can safely be skipped by the reader familiar with this material. 

We first describe the main concepts which will figure prominently in the estimates 
to follow. After this, we proceed to a more detailed review of the features which will 
be especially important for analyzing the cosmology of F-theory GUTs. To this end, 
in the next subsection we provide a brief introduction to the FRW Universe, next 
describing the history of the standard cosmology as the Universe evolves in such 
a Universe. As will be apparent in later sections, the cosmology of the gravitino 
plays an especially important role in assessing the viability of a supersymmetric 
model. For this reason, we next review the computation of relic abundances, and 
in particular, describe in some detail the usual "gravitino problem". While this is 
indeed a significant constraint on many models, in certain situations the dilution due 
to a late decaying cosmological modulus can significantly alter this analysis. With 
this in mind, we next review the primary consequences of late-decaying cosmological 
moduli. We next review the constraint on the number of thermalized relativistic 
species derived from BBN. Finally, we conclude this section with a discussion of 
more refined features of BBN which are present in models such as F-theory GUTs 
which possess a late-decaying NLSP. 
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2.1 FRW Universe 



Below the scale of compactification, and before the present epoch, to good approx- 
imation, the early Universe is described by the four-dimensional Robertson- Walker 
(RW) metric: 

ds^ = dt" - a\t) (y^^ + [de^ + sin^ Odcf")^ , (12) 

where in the above, we have introduced the scale factor a(t), as well as the curvature 
constant fc, which after a suitable rescaling takes the values k = +1,0,-1 for a 
respectively closed, flat, or open Universe. This describes the evolution of an isotropic 
Universe with homogeneous energy distribution. The overall expansion rate of the 
Universe is measured by the Hubble parameter: 

(13) 

It is also common to introduce a dimensionless variant of H, called h defined by the 
equation: 

H = lOO/i km Mpc^^ sec"^ . (14) 
The present value of h is given by [21]: 

ho ~ 0.7. (15) 

When the context is clear, we will often drop this subscript to avoid cluttering various 
equations. 

The background energy density determines the expansion rate of the Universe 
via the Friedmann equations: 

4:71 

d = — —GN{p + Sp)a (16) 

where in the above, p denotes the pressure of the given "fluid", Gn denotes the 
four-dimensional Newton's constant, and p corresponds to the total energy density 
of the Universe. The critical density pc is defined as the value of the total energy 



9 



density for which /c = 0, so that: 



Pc = (18) 

Note that the critical density is a non-trivial function of t. 

The total energy density will in general receive various types of contributions, so 
that p is given by a sum of the form: 



P 



Y.P^- (19) 



Plugging this expression into equation f[T7|) now yields: 



"« = E". = i + 7?l? = ^ + l (20) 



where we have introduced the parameter: 



a = -. (21) 

Pc 

The sign of k correlates with the magnitude of ^tot'- 

l^tot > 1 : A; = +1 (22) 
fitot = 1 : A; = (23) 
fitot < 1 : A; = -1. (24) 

Due to the overall dependence on in pc, it is sometimes convenient to introduce 
the quantity VLih?, with h = kH with k the constant implicitly defined by equation 
fini) . The overall dependence on H"^ factors out, yielding: 

n,h^^P^.{.Hf = ^-^^.p,. (25) 
Pc 3 

The energy density of the i*^ species also evolves with scale. In the approximation 
where the energy density pi is proportional to the pressure pi, the equation of state 
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is given by: 

Pi = WiPi, (26) 
for some constant Wi. The scaling behavior of p for three common choices is: 

Pr oc : Wr = 1/3 (27) 
Pm (X : Wm = (28) 
Pa oc const : w\ = —1 (29) 

where r, m and A respectively denote "radiation" or relativistic matter, non-relativistic 
matter, and vacuum energy density. 

Observation indicates that the Universe has recently transitioned from an era 
of matter domination to one where a background vacuum energy density plays a 
dominant role, with [21]: 

~ 0.7 (30) 
~ 0.3. (31) 

The matter content Qm further subdivides into a subdominant "visible" matter con- 
tribution with ^visible ~ 0.05) with the rest being comprised of the so-called dark 
matter Q^m ^ 0.25, which by definition interacts weakly with the Standard Model 
degrees of freedom. The existence of such a large additional component of matter 
indicates that the Standard Model must be extended in some fashion. An important 
feature of this is that the overall energy density is quite close to the critical value: 

fitot = f^A + ~ 1- (32) 



2.2 Timeline of the Standard Cosmology 

Having reviewed the main concepts which we shall use throughout this paper, we 
now describe in reverse chronological order the timeline for the standard cosmology. 
At various points we also indicate where deviations from this trajectory are possible. 
This material can be found in many standard textbooks, such as [22,23], for example. 
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Era of Structure Formation: T ~ 1 meV^ 100 meV At present, the 
Universe is roughly 10^*^ sec old and the background photon radiation is at a tem- 
perature of 2.7 K ~ 1 meV, which is presently characterized as an era where large 
scale structures have formed. In addition, the energy density of the Universe is 
composed of roughly 5% visible baryonic matter, 25% non-baryonic matter, or "dark 
matter", and 70% of some other type of energy density which appears to share the 
characteristics of a vacuum energy density, or cosmological constant. At this point, 
the abundances of all of the light elements have been produced via cosmological 
processes earlier in the history of the universe, as well as other, more astrophysical 
processes having to do with the birth and death of stars, for example. At somewhat 
earlier timescales around 10^^ sec described by a temperature of a few hundred meV, 
the present galaxies begin to form. 



Era of BBN and Atom Formation: T ~ 1 eV^ 5 MeV Prior to the 
start of structure formation, radiation and matter decouple at roughly 10^^ sec and 
a temperature of ~ 1 eV. At this point, all of the present atoms have formed and 
will soon begin to form large scale structure. At only slightly earlier times, or 
higher temperatures, the universe is sufficiently hot to overcome the binding energy 
of atoms. At this point, nucleosynthesis has finished, and an appropriate abundance 
of the ions such as H^, D^, T"*", ^i/e"*"^, ^He~^^, '^Li~^^~^ have formedE These 
ions will soon combine with electrons to form neutral atoms. Heavier elements are 
produced later as a result of astrophysical phenomena. Other than '^Li, the observed 
abundances of these elements agree quite well with theoretical computations. As we 
shall later review, this information is sufficiently accurate to provide a strong bound 
on the number of additional relativistic species in thermal equilibrium. BBN begins at 
a temperature T ~ 1 — 5 MeV. Generating the appropriate light element abundances 
requires a sufficient number of baryons. Letting ub and denote the respective 
number densities of baryons and anti-baryons in a comoving volume, BBN requires 
that the ratio between the net baryon number density and photon number density, 



is: 



■^As reviewed in [22], for example, astrophysical processes also account for a fraction of the '^He 
abundance. 
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Era of Radiation/Coherent Oscillation Domination: T ~ 1 GeV^ Tb 

Before the start of BBN, the temperature of the universe is sufficiently hot that 
the particles of the Standard Model exist in a plasma. In standard cosmological 
scenarios, it is assumed that some initial amount of baryon asymmetry necessary for 
BBN has been generated at a temperature Tg. At this point, the particles generated 
during baryogenesis are in thermal equilibrium in a hot plasma. As the Universe 
cools and expands, a particle species may fall out of thermal equilibrium. 

The abundance of these species is particularly important for cosmology. The 
essential point is that each frozen out species will contribute to the present matter 
abundance Om^^- In particular, because the observed dark matter abundance is 
given by QoMh"^ ~ 0.1, each species must satisfy the bound: 

Qih'^ < 0.1. (34) 

When this inequality is not satisfied, it is common to say that the species "overdoses" 
the Universe. As explained in [22], the more precise statement is that when Qih'^ is 
too large, the predicted value of h would not be in accord with observation. 

In extensions of the standard cosmology, additional dynamics beyond those as- 
sociated with the freeze out of various species may also play an important role. For 
example, at temperatures above the start of BBN, scalar fields undergoing coherent 
oscillations can come to dominate the energy density of the Universe. When such 
particles decay, they can release significant amounts of entropy into the Universe. 
This has the effect of diluting the abundance of various frozen out species. In addi- 
tion, the decay products of such fields can also increase the relic abundance of certain 
species. 

Era of Baryon Asymmetry Creation: T ~ Tb — ^> Tp,fj At a temperature 
Tb, the baryon asymmetry requisite for successful BBN is assumed to be created. As 
found by Sakharov, generating an appropriate baryon asymmetry requires that first, 
the accidental U{1)b baryon symmetry of the Standard Model must be broken, and 
further, that both C and CP must be violated. In addition, the baryon asymmetry 
must be generated due to some out of equilibrium process. In order to achieve the 
required value of rjB, there must be an additional source of CP violation beyond 
that which is present in the Standard Model. In scenarios where Tg is close to the 
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GUT scale, off diagonal gauge bosons of the GUT group can potentially generate the 
required values of rjB- Lower temperatures for Tb are also possible depending on 
the particular mechanism in question. For example, in leptogenesis, an asymmetry 
in lepton number is converted to an asymmetry in baryon number via sphaleron 
processes. In these examples, heavy right-handed Majorana neutrinos decay to 
Standard Model particles, generating the initial lepton number asymmetry so that 
the associated temperature Tb is roughly given by the Majorana mass of the right- 
handed neutrinos. An important consequence of this is that the initial reheating 
temperature T^^ must satisfy the inequality T]j^ > min (Mmaj), where Mmaj is 
shorthand for the Majorana masses of the heavy right-handed neutrinos. As we will 
review later, this can lead to a certain amount of tension in many supersymmetric 
models which typically aim to lower to avoid overproduction of gravitinos. Lower 
values for Tb are also possible in less standard leptogenesis scenarios, as well as in 
models where the coherent oscillation of a field generates a large lepton or baryon 
number. This latter scenario, known as the Affleck-Dine scenario is particularly 
attractive in models where other cosmological constraints require Tb to be relatively 
low compared to the value required for other baryon asymmetry scenarios. 

Era of Speculation: T ~ T^^ Mpi Above temperatures where the baryon 
asymmetry is created, there is some initial temperature T^^j corresponding to the 
"reheating" (RH) of the Universe. Below this temperature, radiation domination 
commences. This temperature is typically associated with the end of some era of 
where some high scale dynamics generates the required density perturbations and 
flatness of the present Universe. We stress that below the temperature T^^, the 
primary mechanism which sets these initial conditions, which could originate from 
a mechanism such as inflation or string gas cosmology is inconsequential for the 
analysis of the paper. In this way, in adhering to the philosophy of [1-5] these issues 
can be deferred to a later stage of analysis. 

2.3 Thermodynamics in the Early Universe 

In the previous subsection we provided a rough sketch for the evolution of the stan- 
dard cosmology. We now discuss in greater detail the thermodynamics of the Uni- 
verse, and in particular, review the computation of the relic abundances for hot and 
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cold relics. 



2.3.1 Equilibrium Thermodynamics 

The expanding Universe corresponds to the stage on which the interactions of a given 
particle physics model will play out. We now summarize some basic features of equi- 
librium thermodynamics in the early Universe. Much of the following discussion is 
explained in lucid detail in Chapter 3 of [22]. Our aim here is to give a rough intuitive 
summary of the various formulae which will be important in later discussions. 

Assuming that some high scale dynamics sets the initial temperature of the Uni- 
verse, which we denote by T^h^ can follow the subsequent evolution of a cos- 
mological model. At sufficiently high temperatures, various particle species will be 
in thermal equilibrium. The corresponding interaction rates Tint within the thermal 
bath are specified by the collision time for the particle species in question so that: 

Tint ~ rii {aiVi) (35) 

where (aiVi) denotes the thermally averaged cross section for the species and rij 
denotes its number density. In the limit where the temperature is much greater 
than the chemical potential, the number density, and energy density of a relativistic 
species (m ^ T) are given by: 

rirei ~ (36) 

Prel ~ (37) 

for a relativistic species. In the other limit where the chemical potential /i dominates, 
the above formula holds with T replaced by fi. For a non-relativistic species of mass 
m, the number and energy density are: 

.3/2 / m- ^ 
) exp ( 

pn-re\ ~ mn^.j-eh (39) 

Finally, the entropy density s of the thermal bath is primarily determined by the 



?^n-rel ~ {mT) CXp ( ) , (38) 
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interactions of the relativistic species: and is given by: 



s ~ g.siT)T', (40) 

where here, 

,.,(T)^5^,,(§)%I$:,.(^)' (41) 

bose ^ ' fcmii ^ ' 

with gi the number of internal degrees of freedom associated to a given species, so 
that for example, an electron and positron both have (?e+ = fi'e- = 2. Further, Tj 
denotes the actual temperature of the given species, which in general may differ from 
T. Nevertheless, this distinction is largely unimportant when a given species is in 
equilibrium with the background bath. It is also convenient to introduce a count of 
the total number of relativistic species defined as: 



bose ^ fcrmi 



T 



(42) 



We note that at high temperatures, Tj ~ T, and g^{T) and g^s{T) may be used 
interchangeably. In the high temperature limit where all degrees of freedom are 
relativistic, the value of g^, in the Standard Model and MSSM are respectively: 

g,{SM) = 106.75 (43) 
g,{MSSM) = 228.75. (44) 

In the present era, the net energy and entropy density are given, for example, in 
Appendix A of [22]: 

Pe,o ~ (8.1 X 10"^^) GeV^ (45) 
So ~ 2.3 X 10"^^ GeV^ (46) 

where in the above, the subscript reflects the evaluation of this quantity at present 
times. 

The connection between the expansion of the Universe and the temperature of 
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the thermal bath follows from the fact that the total entropy in a co-moving volume: 



S oc [g^sT^) ■ «^ = constant. (47) 

As a resuh, the temperature of the universe evolves as: 

T^g;s'^\-\ (48) 

Finally, in an era of radiation domination, the fact that the energy density scales 
as p oc 5'*T^, in tandem with the second Friedmann equation fll7p . implies that the 
Hubble parameter is related to the temperature as: 



H'-9.^. (49) 

This relation will be quite important when we discuss the decoupling of thermal 
relics during an era of radiation domination. 



2.3.2 Relics and Decoupling 

In the above, we have implicitly assumed that all species in question remain in 
thermal equilibrium. In this subsection we review what happens when a species 
decouples from this thermal bath. Much of the material of this subsection is reviewed 
in greater detail in Chapter 5 of [22], and we refer the interested reader there for 
further discussion. 

As the Universe expands, the thermal bath cools, and a given species may de- 
couple. This occurs when the associated comoving volume becomes too large to 
allow efficient interactions, so that the i^'^ species "freezes out" at the temperature 
T/ implicitly defined by: 

H{T/) ~ m (aiVi) . (50) 

More precisely, the evolution of the number density Tli clS db function of t is given by 
the Boltzmann equations in the presence of a dissipation term which accounts for 
the overall expansion of the Universe: 

^ + SHm = a, (51) 
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The left hand side corresponds to the time evolution of the number density, and Cj 
is determined by the reaction rates of the thermal bath which can generate the i*^ 
species. The principle of detailed balance implies that Ci is given by: 



Ci = {aiVi) (n^EQ 




(52) 



where nj^EQ denotes the equilibrium number density of the i species. 

After a species is no longer in contact with the thermal bath, its number density 
redshifts with the expansion of the Universe, scaling as a~^. On the other hand, 
returning to equation fHTl) . it also follows that when the entropy remains constant 
within a comoving volume, that the entropy density s will also scale as a~^. In 
determining the relic abundance associated with a frozen out species, it is therefore 
convenient to introduce the yield: 



which modulo subtleties connected to changes in the entropy, remains constant after 
the i*'* species has frozen out. During an era of radiation domination t oc T~^, and 
the Boltzmann equation for the yield attains the form: 



where in the above, we have introduced the parameter: 



After a species has frozen out, the left hand side of equation ( !54l) is to leading order, 
neglig ibleil Hence, the yield at the time of freeze out Yi^oo is given by I^.eq evaluated 
at the freeze out temperature: 



Using the yield, we can determine the relic abundance of the i species. The key 
point is that because the yield does not change after freeze out, the number density 

■^See Chapter 5 of [22] for a more precise discussion based on integrating the Bohzmann equations. 




(53) 




(54) 




(56) 
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at present times is given by: 

ni,o = So ■ Fi,oo- (57) 
As a consequence, the relic abundance is: 

ilih = — = ■ = rriiYi^oo ~ 2.8 x 10 GreV • niiYi^oo, (58) 

Pc,0 Pcfi Pc,0 

where in the final equality, we have used the explicit values of pc,o and Sq given by 
equations (145|) and (146|1 . 

The actual yield of the i^'^ species strongly depends on whether it is relativistic, 
or only semi-relativistic at the time of freeze out. In the former case, the freeze out 
temperature is far above the mass of the given particle, so that it is appropriate 
to use the number density of equation fl36p . In the latter case, the mass may be 
comparable to, or larger than the freeze out temperature, in which case the number 
density is given by equation fl38l) . Restoring all numerical factors, and using the 
value of the entropy density in equation fHU]) . the yield at the time of freeze out for 
a relativistic species is: 

„rel 

yrel ^ i,oo ^ 9cS 

s{Tl) g,s{xi) 

where in the above, Qes = 9 for bosons and g^s = "ig/^ for fermions, with g the 
number of degrees of freedom associated with the given species. 

The evaluation of the yield in the case of a species which is at most semi- 
relativistic at freeze out is somewhat more involved. Nevertheless, for our present 
purposes, the main point is that in a rough approximation, the number density is 
given by evaluating the equilibrium number density at the temperature of decou- 
pling. Returning to the freeze out condition of equation fl50|) . the number density at 
freeze out is: 

<^ - ^ - (60) 



Dividing by the entropy density at the time of freeze out, the yield is therefore: 

n-rel 1/2 / 

yn-rcl _ "'i,oo 1 g,, {1^ ) 

~ s{T/) ^ MpL g,s{T/) m^ {^^V^)' 

Plugging back into equation ( l58i) . and restoring all numerical factors, the relic abun- 



19 



dance is then given by: 



nfh^ ~ 8 X 10 



-2 _ 5'cff (rriiX 




(62) 



10^- 



MpL g,six{) {(TiViY 



(63) 



As a point of terminology, relics which decouple when they are relativistic are often 
called "hot" , whereas relics which decouple when they are non-relativistic are called 
"cold" . An important feature of the cold relic density is that it is inversely propor- 
tional to the thermally averaged cross section {(TiVi). Assuming that (o"jWj) ~ a'^/M'^ 
for a a fine structure constant on the order of ~ 1/50, in a model such as the MSSM 
or Standard Model where ~ 100, a suggestive feature of the above formula is that 
nfh'^ ~ 0.1 when: 



An important caveat to the above computations is that it implicitly assumes that 
the Universe starts at a high enough temperature that the species in question is in 
thermal equilibrium, and then falls out of equilibrium. Indeed, it is in principle also 
possible to consider scenarios where the production of a given species is truncated 
because of T]j^ being lower than the freeze out temperature of the species. Although 
somewhat ad hoc, this is one mechanism which has often been invoked to avoid over- 
production of gravitinos in models where this particle is the LSP. 

2.4 The Gravitino and its Consequences 

In the context of F-theory GUTs, the gravitino corresponds to the LSP. Due to the 
fact that R-parity is typically preserved in such models, the gravitino is stable, and 
can potentially correspond to a cosmological relic. For example, in the context of 
high scale gauge mediation scenarios, the gravitino can have a mass as high as 1 
GeV, although in the specific context of F-theory GUTs, this value is somewhat 
lower at 10 — 100 MeV. In this subsection we review the fact that in many models, 
the gravitino relic abundance can overdose the Universe. Indeed, one of the aims of 
the present paper is to explain how F-theory GUTs naturally solve this "problem" . 



M 



1 TeV. 



(64) 
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2.4.1 Freeze Out of the Gravitino 



Because it interacts so weakly with other particles, the freeze out temperature Ty^ 
of the gravitino is typically quite high. To estimate the value of Ty^y we first clarify 
how the gravitino interacts with the thermal bath of MSSM particles. Following 
the discussion in for example [24], after supersymmetry is broken, the longitudinal 
mode of the gravitino xp^^ eats the spin 1/2 Goldstino mode, -j/^ associated with 
supersymmetry breaking so that: 

—d"^. (65) 

Labeling the bosonic component first, given a chiral multiplet x) vector mul- 
tiplet (A^, A) the Goldstino mode couples to these fields through the associated 
supercurrent so that the gravitino Lagrangian density contains the terms [25]: 

L3/2 D ^ [7^ 7l i^^F^u + ^ ^\<P* + h.c, (66) 

m3/2 • MpL ms/2 ■ MpL 

where in the above, the m's denote the masses of various particles, the 7's are the 
usual Dirac matrices, and we have dropped various constants which are not crucial 
for the discussion to follow. Letting m denote the characteristic mass scale associated 
with the mass splitting between members of a given supermultiplet, it follows that 
the relevant cross section for the gravitino is of the form: 

^3/2 ~ -rj2- {-^] . (67) 

When in equilibrium, the primary thermal production mechanism for gravitinos 
is given by the conversion of particles of supersymmetric QCD into gravitinos via 
processes of the form: 

AB ^ C^3/2, (68) 
where here. A, B, C are shorthand for quarks, squarks, gluinos and gluons so that: 

^3/2~^f^y, (69) 
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where rrig is the mass of the gluino. We refer the interested reader to [24,26,27] for a 
complete hst of interactions and the detailed form of the corresponding amplitudes. 

Returning to cosmological considerations, the gravitino roughly freezes out at the 
temperature Ty^ defined by: 

H{Ti/2) ~ ^3/2 (t^3/2f3/2> • (70) 

Precisely because the gravitino only interacts quite weakly with the background 
thermal bath, it decouples when it is still relativistic. For this reason, it is appropriate 
to use the relation 713/2 ~ for a relativistic species. Furthermore, because the 
decoupling of the gravitinos happens at high temperatures, this decoupling occurs 
in an era when radiation dominates the energy density of the Universe. We note 
in passing that after the gravitino decouples, there could be a transition to a more 
exotic epoch where matter, or the coherent oscillation of a field dominates the energy 
density of the Universe. Using the relation between temperature and the Hubble 
parameter in an era of radiation domination provided by equation poj) : 



(71) 



with the total number of relativistic degrees of freedom, it follows that the freeze 
out temperature satisfies: 

Tl.-9l^'M,,(^]\ (72) 



nig 



where in the above, we have set V3/2 ~ Ij as appropriate for a relativistic species. 
Including all appropriate numerical factors and performing a more precise estimate 
based on integrating the Boltzmann equations yields a value Ty^ ~ 10^° — 10^^ GeV 
for a gravitino of mass 7713/2 ~ 10 MeV [28,29]. Indeed, in comparing the overall 
gravitino relic abundance obtained for T^^^ < Ty^ (a computation we will shortly 
review) with the value in the opposite regime where T^j^ > Ty^, continuity of the 
gravitino relic abundance across this interpolation yields: 

T/, ^ 2 X 10- GeV • (^^Y ( ^ V • (73) 
3/2 V 10 MeV/ ^ - I ^ ^ 



nig 
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2.4.2 Gravitino Relic Abundance 



Having specified the temperature at wliicli tlie gravitino freezes out, in tliis subsection 
we determine tlie corresponding relic abundance. As in previous subsections, at this 
point we will not assume that any late decaying field dilutes the total entropy of the 
Universe. Assuming that the initial reheating temperature T'^fj is greater than the 
freeze out temperature, the formulae for the relic abundance of a "hot" species given 
by equation (l62l) is applicable so that: 

> 4, : nl„k- ~ 8 X 10- . -jy (!^) , (T4) 

where in the above, denotes the fact that these gravitinos are thermally pro- 
duced. Since the gravitino decouples at such high temperatures, the value of g^si^i) 
is given by the total number of degrees of freedom in the MSSM. Including all relevant 
numerical factors, the resulting relic abundance of gravitinos is [30]: 

A perhaps distressing feature of this formula is that for a gravitino of mass m3/2 > 0.2 
keV, the relic abundance would appear to overdose the Universe! 

As the above analysis shows, if the MSSM thermal bath starts out at a very high 
temperature and only cools to the gravitino freeze out temperature at some later 
stage of expansion, there is a risk that the abundance of gravitinos could overdose 
the Universe. For this very reason, it is common in the literature to consider scenarios 
where thermal production of gravitinos has been truncated by lowering the initial 
reheating temperature T^fj below Ty^. 

To determine how low T^^ must be in order to avoid an over-production of grav- 
itinos, we next repeat our analysis of the freeze out temperature of a species detailed 
in section 12.3.21 but now in the more general case where the start of thermal produc- 
tion commences either above or below the freeze out temperature of the gravitinos. 
The main change from our previous analysis is that the relic abundance is now deter- 
mined by the yield at the temperature T^//, rather than the freeze out temperature. 
As throughout this review section, our aim here is to give a rough derivation of this 
formula. More precise derivations may be found for example in [24,26,27]. 
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To estimate the abundance of thermally produced gravitinos, we again use the 
fact that = n^j2/s is roughly constant after the initial production of gravitinos. 
Here, it is important to note that in principle, the initial reheating temperature Tj^^j 
can either be greater than, or less than the temperature at which gravitinos freeze 
out. In the latter case, the thermal bath of MSSM particles will begin producing 
gravitinos up until the temperature at which they freeze out. On the other hand, 
in the latter scenario, the scattering processes described above will convert MSSM 
particles into gravitinos at a temperature, T*^^ which immediately freeze out. In this 
case, the yield of gravitinos is given as: 

Yy, {T,,) - ^ ^ i7(TOJs(ro,)' ^''^ 

where the final estimate follows from the Boltzmann equation ( l5Ti) . In this case, the 
equilibrium number density appearing in C3/2 as in equation is given by the number 
density of the background radiation ririT^f^) so that: 

C3/2(T°^)~(a3/2^^3/2>n^(T°^). (77) 

Using the relation ririT^jj) = s(T]j^), the yield is: 

n% inn) - 91'\TIh) (cT3/2t'3/2> MpA. (78) 
Utilizing our expression for the cross section in equation (!69l) . we find: 

As a consequence, the relic abundance when T^jj < T.{^,^ is: 



-3/2 

1/2 (rpO \ ^RH '"'g 



MpL 1713/2 



^0) 



Restoring all numerical factors as in [27]: 



RH - -3/2 • "3/2- - - ^^Qio GeV; V ras/2 J 



/ VI TeV/ 
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This is to be contrasted with the rehc abundance of gravitinos obtained in the case 
where T^^ > where essentially the same formula applies with T^^^ replaced by 
the freeze out temperature Ty2- Defining: 

T37- = min(T°^,T//2), (82) 

the gravitino relic abundance in either case is given by: 

2 7 x 103 ( \ ( 10 MeV \ / m, y 

Note that when > T^^^, substituting T^^^ of equation fl72p into the above relation 
reproduces the earlier form of the gravitino relic abundance given by equation (174|) : 

nn > Til, : ^,i,h- ~ (84) 

Equation ([HSD illustrates the general puzzling feature of models with a light, stable 
gravitino. On the one hand, particle physics considerations have a priori nothing at 
all to do with the value of T^^f. On the other, the overclosure constraint: 

^3/2/1^ < 0.1 (85) 

is badly violated for m3/2 > 2 keV when > T^j,, and when T^^^ < T^^, imposes 
the stringent condition: 

when the gluino has mass m-g ~ 1 TeV. 

One of the remarkable features which we shall find in section H] is that these sharp 
bounds are significantly weaker in F-theory GUT models. In fact, taking the most 
natural range of parameters dictated from purely particle physics considerations, we 
find that this constraint is completely absent in F-theory GUTsl This is due to the 
fact that the decay of the saxion can release a significant amount of entropy, diluting 
the relic abundance of the gravitino. 
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2.5 Cosmological Moduli and Their Consequences 

In the previous subsection we alluded to the crucial role which the decay of the 
saxion can play in F-theory GUT scenarios. Here, we review the main effects of late- 
decaying moduli for cosmology. To begin, we clarify our nomenclature for "moduli" . 
In the cosmology literature, it is common to refer to any scalar field which has a 
nearly flat effective potential as a "modulus". Prominent examples in the context 
of supersymmetric models are field directions which are massless in the limit where 
supersymmetry is restored. In the context of string constructions, however, moduli 
fields typically refer to deformations of the metric or a given vector bundle which 
can in principle be stabilized by high scale supersymmetric dynamics. To distinguish 
these two notions of "moduli" we shall always refer to fields which develop a mass 
due to supersymmetry breaking effects as "cosmological moduli" . In this section we 
show that the coherent oscillation and subsequent decay of cosmological moduli can 
have important consequences for cosmology. 

Cosmological moduli are commonly thought to pose significant problems for cos- 
mology. Indeed, as we now review, such late decaying fields if present can come to 
dominate the energy density of the Universe. If such fields decay too late, they will 
disrupt BBN. On the other hand, there is also the potential for such moduli to re- 
solve various problems via their decays. This fact in particular will prove important 
when we turn to the cosmology of F-theory GUTs. 

By definition, the effective potential for cosmological moduli are quite flat. After 
the initial reheating of the Universe ends at a temperature T^/^, the correspond- 
ing scalar modes will have a non-zero amplitude, which we denote by (po- As the 
Universe cools, these cosmological moduli develop an effective potential, and begin 
to oscillate about their minima. In particular, whereas in the standard cosmology, 
radiation dominates in the range of temperatures between T^^j and the start of BBN 
Tbbn ~ 10 MeV, when the initial amplitude of a modulus field is large enough, 
this contribution can come to dominate the energy density of the Universe. We now 
review the estimate for determining when this can occur, and also elaborate on some 
consequences associated with the subsequent decay of such a modulus field. 

Letting V{(j)) denote the effective potential for the modulus, the equation of 
motion for the field is given by: 

^ + 3H^ + V'{(f)) =0 (87) 
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where the dots above denote derivatives with respect to time and the prime denotes 
the derivative with respect to cf). For simphcity, we shall consider the special case 
where ^(0) = m00^/2. In this case, the energy density stored in the field is given 
by: 

= + \ri\<^\ (88) 

where in the above denotes the mass of the modulus in question. In the following 
analysis we shall always assume that the mass is roughly constant as a function 
of temperature. When we discuss the cosmology of oscillating axions, we will see 
that temperature effects play a more significant role. 

There are in principle two possibilities for when a modulus will begin to oscillate. 
The first possibility is that below the initial reheating temperature, the effective 
potential may develop below the initial temperature of reheating, in which case the 
modulus begins to oscillate during an era of radiation domination. On the other 
hand, it is also possible that the modulus could begin oscillating during an era above 
that set by the initial reheating temperature. 

Assuming that the modulus begins to oscillate during an era of radiation domi- 
nation, the temperature at which the Hubble parameter becomes comparable to the 
mass scale determines the oscillation temperature so that: 

iJ(rt) ~ m^. (89) 

Using the radiation domination relation H ~ gV^T'^ / MpL, the resulting oscillation 
temperature is: 

Restoring all numerical factors, a more exact analysis yields: 

Tt~0.3v/^V^^- (91) 

We shall frequently refer to this temperature as the "oscillation temperature" of the 
modulus field. 

At temperatures where the modulus field begins to oscillate, the energy density 
p<^ stored in the modulus field is governed by the initial amplitude of the modulus 
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and the mass of the modulus field so that: 



(92) 



On the other hand, the energy density stored in the background radiation ~ T^, 
so that at the era when the modulus field begins to oscillate, we have: 



Combining this relation with equation (pOj) . it follows that the ratio of and p^ at 
the temperature is: 



An interesting feature of this formula is that 0o; rather than appears. 

Following for example [31], we now determine the conditions required for a modu- 
lus field to dominate the energy density. For simplicity, we consider a scenario where 
a single modulus field undergoes coherent oscillation. This scenario can be general- 
ized to situations where multiple fields oscillate and can all contribute a substantial 
portion of the overall energy density of the Universe. 

Assuming that the Universe is in an era of radiation domination, the transition 
to a modulus dominated energy density can occur provided the energy density stored 
in the modulus, p^ becomes comparable to the background radiation p,,. In other 
words, an era of modulus domination commences at a temperature where: 



Assuming that the oscillation of the modulus eventually dominates the energy density 
at a temperature T^^, the scaling behavior a and pr oc T'^ translates into the 
condition: 




(93) 




(94) 




(95) 




(96) 
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Combined with equation (l94ll . the temperature is given by: 



Tim - -^Tt - vW^^- (97) 

PL PL 



More generally, in scenarios where is either greater than or less than T]j^, 
is given by: 



^dom~lJ-min(Tl,T°^). (98) 

In the above estimate, we have implicitly assumed that the modulus field is suffi- 
ciently long lived that it can dominate the energy density. It is also possible, however, 
that the modulus field may decay before, or after an era of modulus domination is 
able to commence. For example, when (po/MpL is sufficiently low, the resulting value 
of may correspond to a timescale which is longer than the lifetime of (j). In such 
a scenario, the decay of will occur in an era of radiation domination. 

The timescale for the decay of the particle is given by the inverse of F^, the decay 
rate for 0. The modulus will either decay during an era of radiation domination, or 
at the end of domination. In the latter case, we shall assume that the Universe then 
transitions back to an era of radiation domination. In both cases, the corresponding 
temperature T^^^y at the time of decay therefore scales with t^^^^y = F^^ as: 



^decay ^ (^decaj/) ' (99) 

which holds when radiation dominates the energy density of the Universe. The 
relation: 

therefore implies: 

TLay-9:'^'V^W^L. (101) 

Restoring all numerical factors, a more exact analysis yields: 

TLay ~ O.Sv/F^M^. (102) 

We note that if F^ is sufficiently small, this can disrupt the standard predictions of 
BBN, significantly altering standard cosmology. As we now explain, in situations 
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where dominates the energy density of the Universe, it is also common to refer to 
this temperature as the "reheating temperature" of 0, and we shall therefore also 
use the notation: 



r|H = ^la,~o.5v/rvi^. (103) 

In order for the Universe to enter an era of modulus domination, the three tem- 
peratures T^cay '^dom '^RH niust Satisfy the system of inequalities: 

TLay < < min(Tl, T°^). (104) 

Returning to equations (197|) and fll02p . this amounts to a condition on the properties 
of the modulus: ^ 

^^^^L<^^/nl^MpL, (105) 

or: 

< JL no6) 

The decay products of the modulus field can also have significant consequences 
for the cosmology of the Universe. Again, the immediate consequence naturally 
separates into scenarios where the modulus decays in an era of radiation domination 
or modulus domination. In the former case, the primary constraint is that the 
resulting decay products must satisfy the bound (which we shall review later) on the 
total number of relativistic species, in order to remain in accord with BBN. While 
a similar constraint also holds in the case where the modulus comes to dominate 
the energy density of the Universe, another significant effect is the release of entropy 
into the Universe as a result of the decay of the modulus. This increase in entropy 
effectively dilutes the total relic abundance of all species produced prior to this decay. 

Letting Vtf^h? denote the relic abundance of a species z prior to the decay of 0, 
the relic abundance after dilution is then given by: 

^after 

where in the above, s denotes the entropy density, and we have introduced the 



30 



"dilution factor" : 

= (108) 

^after 

A similar analysis to that utilized in estimating the ratio of energy densities in 
equation ( l94l) yields the following estimate for the dilution factor, which can be 
found, for example, in [18,22,32,33]: 

D,-^- % /^^ (109) 



^0 min(ro^,Tot) 



In the above, the minimum of Tjj^^ and T^^^ enters because there are in principal two 
possible scenarios, where either the temperature of oscillation is larger, or smaller 
than T]j^. The essential point is that the lower temperature more strongly determines 
the resulting entropy density of the Universe denoted by She fore, so that it is always 
appropriate to take the minimum of these two parameters. The release of this large 
entropy into the Universe effectively "reheats" the Universe. 

An implicit assumption of equation fll09p is that the dilution factor is less than 
one. Indeed, in scenarios where (f) decays during an era of radiation domination, the 
only effect on the history of the Universe may be an increase in the relic abundance of 
a given species. When is formally greater than one, it follows from equation (I109P 
that the temperature is less than Tf^^, violating one of the implicit assumptions 
used to derive our expression for the dilution factor. 



2.5.1 The Saxion as a Cosmological Modulus 

In supersymmetric models which also solve the strong CP problem via an axion, 
the other real bosonic degree of freedom of the corresponding chiral supermultiplet, 
which we shall refer to as the saxion provides an important example of a cosmological 
modulus. The essential point is that in the limit where supersymmetry is restored, 
the axion and saxion have the same mass. Thus, just as for any other cosmological 
modulus, supersymmetry breaking effects provide the dominant contribution to the 
effective potential of the saxion. 

The presence of the saxion is especially significant in models where the scale 
of supersymmetry breaking is lower than in gravity mediation models. On general 
grounds, many of the moduli of the string compactification in such cases will have 
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large masses determined by high scale supersymmetric dynamics. Indeed, in such 
a scenario essentially the only remaining cosmological moduli are those required 
from particle physics considerations, the saxion being the primary example of such a 
cosmological modulus. The dynamics of this field in particular will play an important 
role in the cosmology of F-theory GUT models. 

2.6 Oscillations of the Axion 

In the previous subsection we reviewed the fact that the oscillation of a modulus can 
alter the evolution of the Universe, leading to an era of modulus domination, as well 
as an overall dilution of all relic abundances due to the release of entropy into the 
Universe. In particular, we also reviewed the fact that the saxion, as a component of 
the axion supermultiplet can play the role of such a cosmological modulus. In this 
subsection we review the cosmology associated with the axion, focussing in particular 
on the effects of its oscillation. 

We first begin by reviewing some details of the axion. By definition, the QCD 
axion a couples to the QCD instanton density so that the Lagrangian density for a 
contains the terms: 

La.^on = y {d^^af + ^e^''P^Trsui3)cF,.Fp^ (110) 

where here, F^u denotes the field strength of SU{3)c, and fa denotes the axion 
decay constant. The field a corresponds to the Goldstone mode associated with 
spontaneous breaking of an anomalous global f/(l) symmetry and takes values in 
the interval: 

- n < a < 71. (Ill) 
Constraints from supernova cooling impose a lower bound on fa so that: 

fa > 10' GeV. (112) 

The upper bound on fa is based on cosmological considerations, and is general more 
flexible. One of the purposes of this subsection is to review the derivation of this 
upper bound. 

The effective potential for the axion is generated by QCD instanton effects and 
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can be approximated using the pion Lagrangian, as for example in section 23.6 of [34]: 

Vaxion{a) = ml {1 - cos a) , (113) 

where ~ 130 MeV is the mass of the pion, and /^r ~ 90 MeV is the pion decay 
constant. The mass of the canonically normahzed axion is therefore: 

„„.!!!j^.6xlO-eV.(12;^). (114) 

At temperatures T ^ ^qcd ~ 0.2 GeV, the mass of the axion depends non-trivially 
on T, and is given by: 



m.(T) = m,-(%^V- (115) 



The axion is a very fong hved particle, and can therefore have consequences for 
cosmology. Returning to the pion Lagrangian, it can be shown that the primary 
decay channel of the axion is into two photons. Following section 23.6 of [34], the 
relative decay rates between a ^ 77 and vr'' — > 77 is given by the square of the ratios 
of the two decay constants, fa and f^, multiplied by an overall phase space factor 
proportional to m^/m^. The relative decay rates are therefore: 



T-\ / n \ Z / \ / f 



:ii6) 



Using the lifetime of ti^ given by r^o ~ 8.4 x 10 sec, it follows that the lifetime of 
the axion is: 

r.-lO-^Sec-f^^^ ' ^ 10^9 sec -fi^^^^-^^ (117) 



fa J \ fa J 

which is far greater than the current lifetime of the Universe ~ 10^^ sec. 

Because the axion is quite long lived in comparison to cosmological timescales, it 
can in principle play an important role in cosmology. Even so, due to its small mass 
and the fact that all couplings of the axion to the Standard Model and MSSM degrees 
of freedom are suppressed by powers of l/fa, the total relic abundance of axions 
produced from thermal processes is typically quite small. We refer the interested 
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reader to [22], for example, for further details on such estimates. 

The primary cosmological issue connected to the axion is the fact that at high 
temperatures, the potential for the axion is nearly flat, and the field can easily be 
displaced from its minimum. Much as for cosmological moduli, the oscillation of 
the axion can then have consequences for cosmology, appearing as a zero momentum 
condensate of non-relativistic particles. The corresponding contribution to the overall 
energy density of the Universe can in principle overdose the Universe, or provide a 
significant component of the overall dark matter. To this end, we now review the 
associated relic abundance from coherent oscillation of the axion0 

2.6.1 Axionic Dark Matter 

As indicated previously, the axion is nearly massless at high temperatures. By in- 
spection of equation (162|) . the thermal production of axionic relics is very small owing 
to the small mass of the axion ~ 10~^ eV. For this reason, axionic dark matter 
is only a viable candidate when produced through some non-thermal mechanism. 
Precisely because of its small mass, the axion can be displaced from its minimum to 
a value such that —n < Oq < vr. Throughout our discussion, we shall assume that 
this initial displacement is given by roughly the same value in causally disconnected 
patches of the Universe. We note that in scenarios where T^^^ < fa, whatever mech- 
anism solves the homogeneity problem will also translate into a uniform value for the 
initial amplitude in the entire causal patch of the Universe. At higher temperatures 
where > fa, the axion is still not well-defined, so that once the Universe cools 
to a temperature below /„, the initial amplitude of the axion may be different in 
distinct patches of the Universe. Nevertheless, this simply amounts to replacing the 
uniform value of the initial amplitude by a rough average over various causal patches. 

^Although it is beyond the scope of the present paper to present speculations on the evolution 
of the Universe at temperatures greater than T^jj, in the context of inflationary models where 
T^jj < fa, quantum fluctuations in the oscillation of the axion of the form a = (a) + Sa associated 
with oscillation of the axion can in principle induce density perturbations, leading to small variations 
in the CMBR. These "isocurvaturc perturbations" occur in models where the fluctuation mode exits 
the horizon during the expected dc Sitter phase and remains frozen in until some time after inflation 
ends, at which point the mode re-enters the horizon. In the context of models in which inflation 
occurs at a temperature Tj^^ < fa, this leads to a bound on the reheating temperature of the 

form: T^j^ < 10^^ GeV • {fl^^h^y^^^^ ■ (/a/lO^^ GcV)^/^^. Sec for example [22] for further details 
of isocurvaturc perturbations. 
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Once the Universe cools sufficiently, the axion will begin to oscillate, creating a 
condensate of zero momentum particles. We now proceed to estimate the effective 
number density of this condensate, and compute the associated relic abundance of 
axions. Just as in our review of general cosmological moduli, under the assumption 
that the mass term dominates the effective potential, the equation of motion for the 
axion is given by: 

'd + 3Ha + ml{T)a = (118) 
where here, we have included the explicit T dependence of given by equation 

(insD. 

A priori, the axion may begin to oscillate during an era of either radiation or 
modulus domination. In the latter case, the decay of the modulus can dilute the relic 
abundance of the axion, so that in principle, the value of the decay constant can be 
increased [18]. However, when this is not done, a comparison of the relic abundances 
obtained from oscillation in an era of radiation domination and the "undiluted" relic 
abundance of the axion obtained from an era of modulus domination are numerically 
quite similar. For this reason, the diluted relic abundance is typically negligible when 
the axion starts oscillating before the modulus decays. Since we are interested in the 
case where the undiluted relic abundances are numerically quite similar anyway, for 
our present purposes it is sufficient to review the relic abundance computation in the 
case where the axion begins oscillating during an era of radiation domination. 

The axion begins to oscillate at a temperature T^^^ where the mass ma{T) is 
comparable to the overall Hubble parameter: 

H^ma{T:j. (119) 

Assuming that it begins oscillating during an era of radiation domination so that 
H ~ gV'^T'^/MpL it follows that T^^^ is given by: 



ql^'^A'^nnn \ fa J 



■ OSC 



Ar 



(120) 



^QCD \gj A^^^^ 

Dropping the weak dependence on fa, the temperature at which the axion begins to 
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oscillate is therefore given by: 

T:^, ~ 10 ■ Aqcd ~ 1 GeV, (121) 

which is only a few orders of magnitude away from the start of BBN, with Tbbn ~ 2 
MeV ~ 10-2 ■ Aqcd. 

We now proceed to determine the relic abundance of the axion. As in our anal- 
ysis of the cosmological modulus, the energy density stored in the axion when it 
commences oscillation is given by: 

Pa(CJ ~ ImliT^Jifaaof. (122) 

Treating the field condensate as a collection of non-relativistic particles at zero mo- 
mentum, the initial number density is: 

MT^sc) - - Irn^iTDifaaof - ^^.(/.ao)^ ■ f • (123) 

The yield of axions is therefore: 

osc) \ ^ osc J 

Including all relevant numerical factors, a similar analysis to the one already pre- 
sented yields the final estimate for the axion relic abundance [23]: 

Thus, under circumstances where is roughly an order one number, it follows that 
overclosure constraints impose the condition < 10^^ GeV, so that axions can in 
principle comprise a component of dark matter. Nevertheless, this is quite sensitive to 
the actual value of Qq so that even when ~ 10"^, the corresponding relic abundance 
will be negligible. Further, in the event that the axion begins oscillating during an 
era of modulus domination, the numerical similarity of the two relic abundances for 
fa ~ 10^2 GeV implies that once the effects of dilution are taken into account, in the 
latter case the relic abundance is always negligible. 
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2.7 Constraints from BBN 



While many extensions of the Standard Model come equipped with potential dark 
matter candidates, it is also quite important to check that any such extension does 
not introduce additional elements which conflict with well-established features of 
the standard cosmology. In this regard, the standard cosmology prediction for the 
abundances of the light nuclei , are in excellent agreement with observation, 
and are in reasonable accord with ^He~^^, '^He^^ . The predicted abundance of '^Li 
derived from the standard cosmology appears to reveal a discrepancy between theory 
and observation. Optimistically speaking, this can be viewed as a potential window 
into the physics beyond the Standard Model which could potentially alter some of 
the reaction rates present in standard nucleosynthesis. 

One of the most remarkable features of BBN is that the resulting abundances of 
light elements essentially depend on only the expansion rate of the Universe, and the 
overall baryon asymmetry: 

VB = — ^, 126 

where in the above, ub, n-^ and respectively denote the number density of baryons, 
anti-baryons and photons. Quite remarkably, although the resulting abundances span 
approximately nine orders of magnitude, they are all correctly accounted for when 
the baryon asymmetry falls within the narrow window: 

4.7 X 10^^° <Vb < 6.5 X 10"^°. (127) 



Extensions of the standard cosmology can potentially threaten this result in one 
of two ways. As we will shortly review, BBN imposes signiflcant limits on the overall 
expansion rate, and as such, effectively constrains the total number of relativistic 
species present at the start of BBN. The presence of late-decaying particles can also 
alter the results of BBN by either destroying, or producing too much of a given light 
element. For example, in the context of the MSSM, in scenarios with a bino-like 
NLSP, the decay of this particle into a photon and gravitino can potentially disrupt 
the production of certain elements through the presence of additional background 
photons. In certain cases, however, such decays can in fact improve the agreement 
between theory and observation. An important example of this type is the overall 
abundance of '^Li, which in the standard cosmology turns out to be a factor of 2 — 5 
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larger than is observed. In the remainder of this subsection, we provide additional 
details on these two central constraints. 

2.7.1 BBN Bounds on Relativistic Species 

As reviewed for example in [22], increasing the expansion rate of the Universe leads 
to an increase in the total amount of "^He produced by BBN. Due to the connection 
between the Hubble parameter and the number of relativistic degrees of freedom 
coupled to the thermal bath: 

H=,mI^T\ (128) 

a constraint on the expansion rate translates into a direct bound on the total number 
of relativistic species. Although it is beyond the scope of the present paper to review 
the derivation of how production of ^He translates into a bound on the expansion 
rate of the Universe, the end result of this calculation is that at the start of BBN, 
the total number of relativistic species must satisfy: 

g,{T ~ MeV) < 12.5. (129) 

where g*(T) as a function of T is given by equation ( 142|) . The relativistic degrees of 
freedom of the Standard Model already nearly saturate this upper bound, and are 
given by the photon {g^ = +2), three species of neutrinos {gi, = 6) and the electrons 
and positrons {g^- = ge+ = 2). Because these particles are all in thermal equilibrium, 
equation simplifies to the special case where Tj = T for all species so that: 

gSM^rp ~ 1 MeV) = 2 + - (6 + 2 + 2) = 10.75. (130) 

8 

Combined with the strong upper bound provided by (11291) . the general temperature 
dependence of g^ provided by equation (l42l) leads to the inequality: 

new bose ^ ' new fermi 



T \ 



1 MeV 



< 1.75. 



(131) 
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Note in particular that even one fermionic species in thermal equilibrium already 
completely saturates this upper bound. For example, this bound implies that only 
one additional species of interacting relativistic neutrinos can be included in an 
extension of the Standard Model, and that in this case, absolutely no additional 
degrees of freedom can be added without disrupting BBN! This same condition can 
also be stated as a bound on the total energy density contributed by an additional 
relativistic species: 

(^) < ^, (132) 

V Pr J BBN 43 

where pextra denotes the energy density stored in the extra relativistic species. 

It is important to qualify that the bound on the number of relativistic species 
is most stringent in the case of additional degrees of freedom which directly couple 
to the background thermal bath. For example, this might appear to contradict the 
possibility of Dirac-like neutrinos, because if additional light states happened to be 
in thermal equilibrium at the start of BBN, inequality f ll3ip would be violated. 

Additional relativistic species could be present if they decouple at a sufficiently 
high temperature, which we denote by To- Comparing the scale factor dependence 
in equations fl27|) and fl28|) with the temperature dependence in equations fl37j) and 
(j39|) . it follows that the temperature of a decoupled species i obeys the relation: 

Decoupled and Relativistic: a oc — . (133) 

On the other hand, species which remain coupled to the thermal bath are more 
directly sensitive to changes in the number of relativistic species. Indeed, equation 
(HHIl implies that the overall scaling of the thermal bath evolves as: 

aT^g;^/\T). (134) 

Comparing equations (I133p and (11341) . the resulting temperature ratio entering in- 
equality (113 ip is: 

T, / 10.75 ^ 



T \ (T M ■ (^35) 
T \g^{Tr))J 

In other words, if a species decouples at a sufficiently high temperature, the actual 
contribution to g^{T 1 MeV) will be significantly suppressed. 
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Returning to the case of right-handed neutrinos mentioned previously, we note 
that such particles decouple at a temperature To such that: 



g*{TD) > 106.75, 



(136) 



where this lower bound corresponds to the number of degrees of freedom of the 
Standard Model. As a consequence, we obtain the relation: 



so that inequality fll3ip remains intact. 

2.7.2 BBN and Late Decaying Particles 

Given the fact that even crude bounds from BBN tied to the expansion rate of the 
Universe translate into detailed constraints on the number of relativistic species, it 
is perhaps not surprising that the reaction rates necessary for generating the correct 
abundance of light elements from BBN are also quite sensitive to the presence of late 
decaying particles. On the one hand, this imposes important constraints on potential 
extensions of the Standard Model, because the abundances of the light nuclei H~^, 
D~^, T"*", ^He~^~^, ^He^^ are all in reasonable accord with observation. On the other 
hand, this also provides a window into new physics, because the standard cosmology 
appears to predict an abundance of ^Lz which is too large by a factor of 2 — 5 when 
compared with observation. We refer the interested reader to [35] and references 
therein for a very recent account of the current bounds on various abundances. 

As briefly mentioned above, late decaying particles are possible in certain super- 
symmetric extensions of the Standard Model. For example, in the context of the 
MSSM where the effects of supersymmetry breaking are communicated via gauge 
mediation, the gravitino is the lightest superpartner of the MSSM, and either the 
bino or stau corresponds to the next to lightest superpartner (NLSP). This is the case 
of primary interest for F-theory GUTs, and so in the remainder of this subsection 
we shall therefore restrict attention to this case. 

The decay rate of the NLSP into a gravitino and its Standard Model counterpart is 
determined by the universal coupling of the gravitino to matter provided by equation 




(137) 
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(|66l) . The calculation of the lifetime of the NLSP is reviewed for example, in [36], 
and leads to the well known result: 



where tjiimlsp denotes the mass of the NLSP and k is a model dependent factor 
which is unity for the case of the stau NLSP, and measures the photino content of 
the bino in the case of the bino NLSP. The particular normalization for the two 
masses has been chosen to conform with natural values in the range expected in the 
specific context of a high-scale gauge mediation model, as is the case in the context 
of F-theory GUTs. The lifetime of the NLSP is to be compared with the timescale of 
BBN, which roughly commences at a temperature oi Tbbn ~ 1 MeV, corresponding 
to the timescale t ~ 0.2 s. By inspection of equation fll38p . it follows that in certain 
situations, the NLSP could potentially decay just before the start of, or even during 
BBN! 

At the most conservative level, the usual results of the standard BBN cosmology 
can typically be retained if the NLSP decays prior to the start of BBN. Returning 
to equation (11381) . decreasing 7713/2 or increasing mi^LSP will both decrease the value 
of tmlsp- In particular, the fact that the fifth power of ttij^lsp appears in t^lsp 
implies that even very mild adjustments in this value can significantly decrease the 
lifetime of the NLSP. 

Assuming that the NLSP decays during BBN, its decay products could poten- 
tially jeopardize the production of the light element abundances, or could bring the 
abundance of light elements such as ''Li into better accord with observation. The 
precise effect of the NLSP depends on whether it decays to a photon and gravitino 
as for a bino-like NLSP; or whether it decays to a tau and gravitino, as in the case 
of a stau NLSP. The impact from a late decaying NLSP depends on details of a 
particular model, such as the overall abundance of the NLSP prior to the start of 
BBN. Nevertheless, under reasonable assumptions such that the relative abundance 
of the NLSP to the overall baryon density is not distorted by dilution effects, it is 
possible to estimate the production of the light elements due to BBN. 

In the case of a bino-like NLSP, the analysis reflected in figure 7 of [37] indicates 
that when the gravitino has a mass of 10 — 100 MeV, the main predictions of BBN 
remain intact. At larger values of the gravitino mass in the range of > 1 GeV, only 




2 
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somewhat large values of the bino mass remain in accord with BBN. Similarly, in the 
case of a stau NLSP, figure 12 of [37] illustrates that the same range for the gravitino 
mass remains in accord with BBN with similar constraints for the mass of the stau 
in the range of larger values of the gravitino mass. 

Interestingly, the very recent analysis of [35] also indicates that a gravitino in this 
same mass range appears to also decrease the overall abundance of ^Lz, bringing the 
resulting abundance into better agreement with observation. Although a complete 
analysis of BBN in the context of F-theory GUTs is beyond the scope of the present 
paper, quite auspiciously, the natural range of parameters of F-theory GUTs suggests 
a gravitino mass in the range of 10 — 100 MeV! 

3 F-theory GUTs and the Axion Supermultiplet 

In this section we study the axion supermultiplet in F-theory GUT models. To this 
end, we first briefly review the main features of F-theory GUTs [1-3] . This is followed 
by a discussion of the axion supermultiplet, and in particular, the interactions of the 
saxion with the MSSM. 

In F-theory GUT models, the singularity type of an elliptically fibered Calabi-Yau 
fourfold with section over subspaces of the threefold base of complex codimension 
one, two and three respectively determine the resulting gauge symmetry, matter 
content and interactions terms of the low energy effective theory. Singularity type 
enhancements along complex surfaces are interpreted as seven-branes with ADE 
gauge group. These seven-branes then intersect over Riemann surfaces or "matter 
curves" where the singularity type enhances by at least one rank. Finally, this 
enhancement can increase further at points of the compactification, corresponding 
to the intersection of at least three seven-branes at a single point. 

Although perhaps contrary to previous experience with string compactifications, 
local F-theory GUTs provide a somewhat rigid framework for string based model 
building. In these models, the existence of a limit where gravity can in principle 
decouple requires that the GUT group seven-brane must wrap a del Pezzo surface. 
Moreover, breaking the GUT group to the Standard Model gauge group is remarkably 
constrained in such models [2,9]. 

Next consider the supersymmetry breaking sector of F-theory GUTs. For the 
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purposes of this paper we shall assume that the primary features of the deformation 
of a minimal gauge mediation scenario developed in [3] are satisfied. The existence 
of a gauge mediation scenario assumes that most moduli present are stabilized due 
to high scale supersymmetric dynamics. Much as in [3], our attitude will be that 
phenomenological constraints on the particle physics content of this class of models 
should be viewed as imposing interesting restrictions on possible global completions 
which satisfy these conditions. 

We now review some further details of the supersymmetry breaking sector dis- 
cussed in [3]. Consistent electroweak symmetry breaking requires that the parameter 
fj, of the MSSM superpotential: 



must not be significantly different from the weak scale. If supersymmetry break- 
ing indeed stabilizes the hierarchy between the weak scale and the GUT scale, this 
naturally suggests that the value of fi should be correlated with the scale of super- 
symmetry breaking. The effects of supersymmetry breaking can be parameterized in 
terms of the vev of a GUT group singlet chiral superfield X such that: 



In [3] , some explicit solutions to the fi problem were obtained under the assumption 
that X localizes on a matter curve which intersects the GUT seven- brane at a point. 
Integrating out the Kaluza-Klein modes on the X field curve generates a higher 
dimension operator in the effective theory of the form: 



where, as estimated in [2], Mx ~ ]^q15.5±o.5 jg ^j^g Kaluza-Klein scale associated 
with the curve supporting the X field and 7 ~ 0(10). In order to generate the 
correct value of the fi term, F must attain the value ~ xo^^^°'^ GeV^. As explained 
in [3], this value is naturally attained via instanton effects associated with Euclidean 
three-branes wrapping other complex surfaces of the geometry. 

Following the model of [3], the supersymmetry breaking sector is closely tied to 




(139) 



{x) = x + e^F. 



(140) 
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the seven-brane associated with the U{1)pq gauge symmetry. This gauge theory is 
anomalous, and as such, instanton contributions can generate contributions to the 
superpotential of the form: 

Winst ^MlQ-q-X (142) 

where: 

g~e-^°'^«. (143) 

Assuming a fixed value for g, this determines the F-term component of X. As 
explained in [3], to properly analyze the PQ symmetry breaking sector, it is necessary 
to treat both q and X as dynamical fields. In that context, it was shown that an 
appropriate tuning in the Kahler potential for q and the fiux-induced FI parameter 
for the PQ gauge theory is compatible with stabilizing the vev of X at the scale 
10^^ GeV. The Goldstone mode associated with the breaking of the accidental global 
U{l)pQ symmetry is parameterized by the phase of the gauge invariant operator 
q - X. In particular, the axion is therefore primarily given by the phase of X, with a 
small contribution from q. The corresponding axion decay constant is then given as: 

fa = V2x ~ 10^2 GeV. (144) 



To a large extent, the allowed mediation mechanism which communicates the 
effects of supersymmetry breaking to the visible sector is dictated by the fact that 
F/MpL is far below the weak scale so that Planck suppressed operators, and therefore 
gravity mediated supersymmetry breaking cannot generate viable soft mass terms. 
By contrast, geometric realizations of minimal gauge mediation scenarios are far 
more viable in this scenario, and explicit models based on minimal gauge mediation 
have been discussed in [3] (see also [11]). In minimal gauge mediation (MGMSB), 
the soft mass terms are completely fixed by the gauge couplings of the MSSM and 
the ratio Fx /x = A as: 

a F a . 

up to numerical factors associated with the representation content of a given field. 
Here, a is shorthand for the fine structure constants of the various gauge groups of 
the Standard Model under which a given superfield may be charged. 

As explained in [3], the F-theory GUT actually corresponds to a deformation of 
minimal gauge mediation. Indeed, precisely because the X field localizes on a mat- 
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ter curve, it will be charged under additional seven-branes of the compactification. 
These seven-branes endow the low energy effective theory with additional, generically 
anomalous U{1) gauge group factors. The generalized Green-Schwarz mechanism 
cancels this anomaly, but the presence of the requisite coupling of the gauge field to 
an axion-like field in the four- dimensional effective theory generates a large mass for 
the gauge boson via the Stiickelberg mechanism. Below the mass scale of this gauge 
boson, the theory will therefore retain an anomalous global U{1) symmetry, which 
in appropriate circumstances can be identified with a U{1) Peccei-Quinn symmetry 
which we denote as U{1)pq. 

The fields of the MSSM are charged under U{1)pq with charges —2, +1 and —4 
for the respective Higgs fields, chiral matter and X field of the F-theory GUT model. 
Integrating out the heavy U{1)pq gauge fields also generates higher dimension oper- 
ators in the low energy effective theory of the form [38]: 

L^ff D -^^exe^ ■ [ d'OX^X^^^ (146) 

where Mu(i)pQ denotes the mass of the heavy U (1) pq gauge boson and the e's denote 
the respective charges of X and MSSM field $ under U{1)pq, and apq is the fine 
structure constant of this gauge theory. Once X develops a vev as in equation 
fll40p . this generates an additional contribution to the masses squared of the scalar 
component of $ at the messenger scale: 

^l,so/t = '^%,MGMSB H i p ■ ^PQ (147) 

\ex I 

where we have introduced the PQ deformation parameter: 



Mu{i)pQ 



(148) 



In terms of the anomalous U{\)pq gauge theory, this contribution can be interpreted 
as a supersymmetry breaking D-term. Insofar as the PQ deformation corrects the 
soft mass terms of the MSSM, the value of ApQ from prior considerations is on the 
order of ~ 100 GeV. As we show later, this estimate is borne out by the cosmology 
of the F-theory GUT scenario. 
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By inspection of equation ( 1147^ . the PQ deformation decreases the soft masses 
squared when e$ - ex < 0. As a consequence, there is a hmit to the size of Apg before 
the PQ deformation induces a tachyon in the squark/slepton sector of the MSSM. 
The precise value of this vev depends on the value of A. For example, in a model 
with a single vector-like pair of messenger fields in the 5 © 5 of SU{5), the minimal 
value of A consistent with the Higgs mass bound rriho > 114.5 GeV is A ~ 1.3 x 10^ 
GeV, and the maximal value of Apq allowed before the mass squared of the lightest 
stau becomes tachyonic is Apg ~ 290 GeV. For vanishing PQ deformation, a bino- 
like lightest neutralino is the NLSP. On the other hand, for large PQ deformation, 
the NLSP can instead correspond to the lightest stau. The PQ deformation also 
plays a significant role in the dynamics of the X field which we analyze in detail in 
subsection 13. 2[ This is particularly relevant for cosmological considerations because 
as we explain in subsection 13. H the axion and gravitino are both closely tied to the 
dynamics of X which is in turn controlled by the value of Apg. 

3.1 Axion Supermultiplet 

In the axion solution to the strong CP problem, an anomalous global U{1)pq symme- 
try is spontaneously broken at an energy scale fa- The associated Goldstone mode 
then corresponds to the axion field, which we denote by a. In a supersymmetric 
theory, a fits into a complete supermultiplet given by one additional real bosonic 
degree of freedom s, a fermionic component ip±, and an auxiliary field F±, which we 
assemble into the chiral superfield: 



The field s corresponds to the "saxion" and the field ip± corresponds to the "axino" . 
By definition, A couples to the QCD superfield strength through the coupling: 



In this section we review and slightly extend the analysis of [3] by showing that the 
components of A are to leading order given by the components of the chiral superfield 
X. In addition to its role in PQ symmetry breaking, the X field also plays a key role in 
supersymmetry breaking. In the field theory limit, this leads to an exactly massless 
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Goldstino mode in the low energy theory. Precisely because X is the primary 
source of supersymmetry breaking, a linear combination given predominantly by 
the fermionic component of the X superfield with smaller contributions from other 
fermionic modes corresponds to the Goldstino. Away from the strict field theory 
limit, the Goldstino is eaten by the gravitino via the super-Higgs mechanism. To 
leading order, the axino therefore corresponds to the longitudinal components of the 
gravitino. 

3.1.1 U{l)pQ Goldstone Mode 

We now describe the Goldstone mode associated with the breaking of the anomalous 
U{l)pQ symmetry. As explained in [3], to leading order, argX corresponds to 
the axion field. Strictly speaking, however, this identification is not completely 
accurate because the axion superfield corresponds to a linear combination of X with 
subleading contributions from other chiral multiplets. The reason for this can be 
traced back to the way in which the PQ gauge boson develops a mass. 

In general terms, the anomalous U{1)pq gauge theory consists of n chiral super- 
fields Xi with charges such that: 



The anomaly for the corresponding U{1) gauge boson is then canceled via the Green- 
Schwarz mechanism. This corresponds to introducing an additional axion-like super- 
field C such that the chiral superfield e'' has charge —Q. C couples to the PQ vector 
multiplet via the D-term: 



where here, K denotes an appropriate Kahler potential. The corresponding gauge 
field develops a mass via the Stiickelberg mechanism, leaving a nearly exact global 
symmetry in the low energy effective theory. Once some combination of the Xj's 
develop a vev, this global U{1)pq will be broken, leaving behind a Goldstone mode. 

The axion supermultiplet is given by a linear combination of C and the associated 
"phases" of the Xj's. Another linear combination of these fields is eaten by the vector 
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multiplet. Assuming a canonical normalization for all of the Xj's, the direction in 
field space fixed by the D-term potential of the PQ seven-brane gauge theory is: 

n 

~ QK' = ^flu., (153) 

i=l 

where here C,fiux denotes a flux induced FI parameter. The directions unfixed by 
the D-term potential are conveniently parameterized in terms of chiral superfields 0^ 
defined as: 

Xj = \xj\exp{iQj) (154) 

The linear combinations of Qj neutral under U{1)pq correspond to possible flat 
directions not fixed by the D-term potential. For example, in the context of the 
Fayet-Polonyi model, the U{1)pq invariant combination: 

X = q-X (155) 

develops a supersymmetry and global PQ symmetry breaking vev. This vev is stabi- 
lized by contributions to the Kahler potential for X and q [3] . The bosonic component 
of G is given by: 

Q = a + is + ... (156) 

where a is the axion and s is the saxion. The mode a takes values in the interval 

—IT < a < IT. 

3.1.2 Supersymmetry Breaking and the Goldstino 

In the above analysis, we assumed that supersymmetry was unbroken. At energy 
scales below x ~ 10^^ GeV, supersymmetry is broken via the Fayet-Polonyi model 
described in [3]. In the strict field theory limit, spontaneous supersymmetry breaking 
generates a massless fermion corresponding to the Goldstino mode. The explicit 
form of the Goldstino mode depends on the details of F- and D-term breaking of a 
particular model. Nevertheless, certain features of the Goldstino mode and how it 
couples to fields of a particular model are universal, and our analysis will for the most 
part only rely on these well known features. As reviewed, for example in [36], the 
explicit form of the Goldstino mode is given as a linear combination of the fermions 
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Xa and Xi respectively from the vector and chiral multiplets: 

G = '-^K + (157) 

where the effects of supersymmetry breaking are encoded in non-zero vevs for some 
subset of the D"^ and F*. Taking into account the fact that Mpi is not infinite, 
spontaneous supersymmetry breaking imphes that the Goldstino mode is eaten by 
the gravitino. The mass of the gravitino is: 

In the exphcit Fayet-Polonyi model presented in [3], there will in general be 
contributions to both the F- and D-term components of the gravitino. In addition to 
contributions to the D-term potential, the Fayet-Polonyi model contains an instanton 
induced contribution to the superpotential of the form: 




W = MIqC^X (159) 

where the mass scale Mpq ~ Mqut which as explained in [3] can in general be 
different from the mass of the PQ gauge boson. For this range of mass scales, it 
turns out that the vev of e'' can naturally attain the value ~ 10~^^, as required to 
achieve F ~ 10^'^ GeV^. With the identification oc we can alternatively 

view this Fayet-Polonyi model as a Fayet model of supersymmetry breaking with 
superpotential: 

W = mXiXn+i. (160) 

As explained in [3], the full sector requires a non-trivial Kahler potential for Xi and 
Xn+i. The F-term components of the various superfields are therefore determined 
as: 

Fi = -mxn+i ~ 10^^ GeV^ (161) 
Fn+i = -mxi = -mxn+i^^ ~ 10^^ GeV^ (162) 

where in the above we have plugged in rough representative values of ~ Mu{i)pQ ~ 
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Mgut and xi ~ 10^^ GeV consistent with the estimates obtained in [3]. 

In addition to these F-term contributions to supersymmetry breaking, as ex- 
plained in [38], we should also expect a contribution from D-term breaking, which 
in the present class of models is given as: 

which is far smaller than the F-term breaking components. 

To summarize the discussion presented above, the gravitino is predominantly 
given by the fermionic component of the X field. Nevertheless, for certain purposes, 
it is important that the gravitino corresponds to a linear combination of the Goldstino 
which contains additional contributions from the fermionic components of e*^. For 
example, in the low energy effective theory for the X field, integrating out the heavy 
U{l)pQ gauge boson generates the higher dimension operator: 

3 / d'eX^XX^X. (164) 

Once X and Fx develop non-zero values, this would appear to induce a mass term 
for the gravitino of order |xApQ/Fjjc|. In a more complete analysis, the fermionic 
mode ipx mixes with the fermionic components associated with . That this must 
be the case follows from the fact that in the limit MpL — » oo, the gravitino is exactly 
massless. Indeed, a naive analysis suggests that the fermionic component ipc has a 
mass term of order |x„+iApq/Fx| which is significantly larger. For most purposes, 
however, this distinction will not play any crucial role in many of the cosmological 
considerations discussed here. 

3.2 Axion Supermultiplet Interaction Terms 

For the purposes of cosmological considerations, it is also important to determine 
the couplings between the axion supermultiplet and the matter content of the rest of 
the F-theory GUT model. Indeed, these interaction terms determine the lifetime of 
the saxion, which can have important consequences if it is sufficiently long-lived. In 
particular, the decay of the saxion can reheat the universe. The specific details of 
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this reheating depends on the dominant mode of decay, as well as the primary decay 
rates. 

As derived in [3j, the effective action for the X field in 4d A/" = 1 superspace 
contains the termso 

Lx = j dWX^ + Rej d-'e^^^^^TrWW^ - j d'OC^ (log \X\'f 

(165) 

- / rfW$XtX - mUX - (166) 

where in the above, $ is shorthand for any MSSM chiral superfield, denotes 
any gauge superfield strength, the multiplicative factors Cw and C$ are determined 
by the gauge couplings constants of the MSSM, as well as the quadratic Casimirs 
for the representations of the messenger fields, and e$ and ex denote the U{1)pq 
charges of $ and X. Here, the specific form of the saxion mass squared term is fixed 
by details of the Kahler potential for X, and the axion-like field C which enters the 
Green-Schwarz mechanism [3]. The explicit coupling of the axion supermultiplet is 
then given by performing the substitution: 

X xe'^ + e'^F (167) 

in the Lagrangian density Lx- Although the PQ deformation dependence of the 
saxion mass is subject to order one tunings depending on the Kahler potentials of 
X and C [3], to make our discussion more concrete, we shall estimate the resulting 
mass using the higher dimension operator: 



- / d'ex^xx^x. (168) 

Ml,,. J 



Substituting in our expression from equation fll67p . we shall take as a rough estimate 



^It is important to note that this effective action is only vahd a temperatures far below the 
messenger scale. Indeed, in the context of models where the dilaton directly couples to the QCD 
field strength, high temperature effects can potentially destabilize the value of the dilaton. Here, 
this problem is avoided because the actual coupling is only present in the low energy effective theory. 
See [39,40] for further details on destabilization of the dilaton in other contexts. Nevertheless, it 
would be interesting to study in greater detail the high temperature phase of this system. 
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the saxion mass squared as: 

ml, = A\ex\AlQ (169) 

or: 

rusax = 4ApQ, (170) 

where in the above we have used the fact that \ex\ =4. 

In the present context, however, the most important interactions are those which 
determine the decay modes of the saxion. The dominant interaction terms originate 
from the self-interactions of the X field chiral supermultiplet, and the soft mass 
terms associated with gauge mediation effects. Expanding the kinetic term |c)^X|^ 
in terms of the axion multiplet yields the model independent coupling between the 
axion and saxion: 

Lx D {fa +fa-s) Q {d,af + ^ {d,sf^ . (171) 
In particular, the decay rate s aa is given by: 

Further details on universal couplings of the axion supermultiplet can be found in [41]. 

Gauge mediation generates soft mass terms for the scalars and gauginos of the 
MSSM, and also induces additional interaction terms for the saxion. The essential 
point is that the soft mass terms for a gaugino. A, and scalar, of the MSSM are 
given as: 

Lsoft 3 -'m\,soft{x)\\ + h.c. - ml,,^f^ {x,x) \(f)f (173) 
where in the above 



F 

mx,soft{x) oc — (174) 



^l,softix,x) oc 



X 
X 



(175) 



As noted, for example, in [32], the explicit x dependence of rrisoft implies that in 
both cases, X, and therefore the saxion directly couples to A and 0. Performing the 
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substitution x t-^ x + X, it now follows that X couples to these fields as: 



IX X 

Lsoft 3 --mx,softix)—\\ + ml^^jt{x,x) — + h.c. (176) 



Assuming that a given decay is kinematically allowed, the decay rate of the saxion 
to MSSM particles is therefore: 



s^MSSM 



2 \ 2 / 2 \ 2 

^soft\ _ 1 f^soft^ 



(177) 



27imsax \ fa J SnApQ \ fa 

which is to be compared with the decay s —>■ aa. The ratio of these decay rates is 



s^MSSM V"^so/t 



(178) 



Although a sharp upper bound is somewhat fiexible, in order to avoid a tachyonic 
mass, it is necessary to assume, for example, that Apg < nisoft ~ 100 — 1000 GeV. 
Introducing the branching ratio to axions 

Bs^aa = (179) 
sax 

where Tgax denotes the total decay rate, we conclude that for typical values of ApQ 
and rusoft, Bs^aa is in the range: 

10-=^ < 5„ < 10-1. (180) 

See figures[T]and[2]for plots of the "toy model" branching ratios Vg^aa/ {NughJ^ s_^mssm+ 
Ts__,aa) as functions of Apg with fa and rrisoft fixed to representative values. Here, 
Niight denotes the number of light MSSM species which are sufficiently light to allow 
a decay channel to an MSSM particle. In realistic models, Nught is a non-trivial 
function of Apg, although for illustrative purposes, we take it fixed in the figures. 

When a decay to an MSSM particle is kinematically allowed, the branching frac- 
tion to MSSM particles will most likely dominate over decays to the axion. This is 
an important constraint in the context of decays to relativistic species, such as the 
axion. In a representative case, we can expect Apg ~ 100 GeV and nisoft ~ 10^'^ 
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Figure 1: Plot of the "toy model" branching ratio of the saxion to axions and one 
(red) to five (black) species of MSSM fields as a function of Apg for fixed values 
of fa = 10^^ GeV and mgoft = 200 GeV. For the decay to a representative MSSM 
field, we have used the crude estimate provided by equation (11771) . This situation 
is somewhat idealized, because as Apg increases, the number of decay channels will 
increase, decreasing the overall branching ratio to axions. 
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Figure 2: Plot of the branching ratio of the saxion to axions and one (red) to five 
(black) species of MSSM fields as a function of Apg for fixed values of fa = 10^^ GeV 
and nisoft = 500 GeV. For the decay to a representative MSSM field, we have used 
the crude estimate provided by equation fll77p . This situation is somewhat idealized, 
because as Apg increases, the number of decay channels will increase, decreasing the 
overall branching ratio to axions. 
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GeV, in which case: 

10-\ (181) 

Restricting to saxion decays to MSSM particles, we note that the primary decay 
channel through MSSM scalars proceed via the Higgs bosons and potentially the 
lightest stau ri for sufficiently large PQ deformation. The saxion can also decay 
to MSSM fermions. Note, however, that in the supersymmetric limit, the axion 
supermultiplet coupling to the MSSM is far weaker. Because of this, our expectation 
is that the saxion couples more strongly to the fermionic superpartners in comparison 
to the fermions of the Standard Model. In the context of the MSSM, the gluinos 
are nearly 1000 GeV in mass. The primary decay channel through MSSM fermions 
therefore proceeds via decays to the two lightest neutralinos Xo ^i-nd Xq, which are 
primarily composed of the bino and wino. The respective masses of these particles 
in the present class of models are m^i ~ 170 GeV and m^g ~ 350 GeV. In all cases, 
the relative branching fractions are roughly determined by the relative masses of two 
species i and j to be: 

r 



*yy 



[182) 



Here it is important to note that the soft mass mgoft niay in general differ from 
the actual mass of the particle. For example, whereas the Higgs field has mass 
rriho ~ 115 GeV, the soft mass term which enters into the saxion decay rate formula 
is more accurately approximated by ~ 600 GeV [3]. 

The saxion also couples to the gravitino. The model independent decay rate 
s V'3/2^/'3/2 is given by: 

1 3 / ^ ^ 



r-^3/.^3/.-^]i|^(^;;^J • (183) 

The relative branching fraction to gravitinos versus axions is therefore given as: 

P _ ^8^-4,^/2^^/2 _ p ^ s^-^3/2'4'3/2 _ 2 f fa \^ f rrisax \ ^ 
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The decay rates obtained above allow us to determine the lifetime of the saxion. 
The lifetime of the saxion is given by the inverse of its decay rate: 



"Tsax — ^ sax ~ ^s^aa ' ^ g^aa^ (186) 

where as before Bg^aa denotes the branching ratio of the saxion to axions. Combined 
with our expression for the decay rate to two axions given by equation fll72p . the 
lifetime of the saxion is roughly given by: 

Tsax = Bg^aaT-giaa ~ B.^aa^ ~ 2 X lO'^ SCC -Bg^aa ("M^^ , (187) 

^PQ V ^PQ J 

which shows that the saxion is long-lived. 

Summarizing, the primary decay channel of the saxion is either to an MSSM field 
such as the Higgs, or to a pair of axions. In terms of the PQ deformation parameter 
ApQ, the relevant decay rates are: 

r„ ~ -% (188) 

Ts-^A/ssM ~ - — 7 ( — 7^ ) (189) 



1 J 




87rApQ V 
1024 A|>Q 


fa ) 
f^PQ 


967r Ml^ ' 





rs-i^3/2V'3/2 ~ 7^772- ( - — ) • (190) 

Having detailed the primary channels of the saxion, we now proceed to the cosmology 
of F-theory GUTs. 



4 Cosmology of F-theory GUTs 

A priori, a seemingly consistent particle physics model could be in severe conflict with 
cosmology. In this regard, cosmological considerations translate into constraints on 
properties of the particle physics model. In this section we study the cosmology of F- 
theory GUTs. Quite remarkably, we find that over the available range of parameters 
dictated by purely particle physics considerations, F-theory GUT scenarios appear 
to non-trivially satisfy many cosmological constraints. Further, we find that much of 



57 



the tension typically present in gravitino cosmology finds a natural resolution in the 
context of F-theory GUTs, which has additional repercussions for issues related to 
generating a baryon asymmetry with a sufficiently high value of Tp,fj. Indeed, many 
popular mechanisms for generating sufficient baryon asymmetry hinge on allowing 
high values for the initial reheating temperature, which is commonly viewed as being 
in conflict with the requirements of gravitino cosmology. This resolution comes about 
because of a remarkable conspiracy in the mass of the gravitino and the expected 
dilution effects associated to the decay of the saxion. Thus, while one component 
of the axion supermultiplet might appear to create a source of tension, the other 
component completely relaxes it away. 

Cosmological constraints also provide an important window into UV sensitive 
features of the F-theory GUT model. Indeed, the PQ deformation is directly sensitive 
to the mass of the anomalous U{1)pq gauge boson. This PQ deformation plays 
a prominent role in the dynamics of the saxion, which can in turn signiflcantly 
impact the evolution of the Universe. In this way, constraints from cosmology on the 
dynamics of the saxion fleld directly translate into seemingly far removed constraints 
on the compactiflcation and particle physics content of the F-theory GUT! 

Because the ensuing discussion has many interlinked parts, we now summarize 
the main features which allow F-theory GUTs to evade the typically stringent bounds 
on the initial reheating temperature derived from bounds on the relic abundance of 
gravitinos. Starting from the initial temperature T^h^ ^^e Universe begins to cool. 
At high temperatures, the associated thermal bath converts MSSM particles into 
gravitinos. Due to the small total cross section of the gravitinos, these particles fall 
out of equilibrium at a relatively high temperature, T.^^^, which as before we will 
denote as the "freeze out" temperature. While all of this is occurring, however, the 
saxion modulus begins to oscillate at a temperature T^^^. The essential point is that 
because its potential is so flat, the saxion will naturally be displaced away from its 
minimum. In fact, for values of the initial amplitude sq of the saxion fleld which 
are quite reasonable from the perspective of F-theory GUTs, the Universe eventually 
enters an era of saxion domination at a temperature TJ^^. At this point, the relic 
abundance of gravitinos is already determined and is given by the estimates already 
presented in section |2l The era of saxion domination continues until the saxion 
decays, at which point it reheats the Universe to a temperature T^^^, releasing a 
large amount of entropy as a consequence of its decay. This dilutes the overall relic 
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abundance of all species, such as gravitinos. 

It remains to say whether the actual relic abundance of gravitinos is sufficiently 
low to remain in accord with overclosure bounds. As we will describe in greater detail 
in the sections to follow, although a priori, the oscillation temperature of the saxion 
and the gravitino mass are unrelated, in the context of F-theory GUTs, typical values 
for the saxion and gravitino mass lead to the relation: 

F-theory GUT : T„^,, ~ 4^,. (191) 

This turns out to have the remarkable consequence that the overall relic abundance 
of gravitinos is independent of Ty^, T^^^ and T^i^fjl Moreover, the relic abundance 
of gravitinos thus obtained does not overdose the Universe, and could potentially 
correspond to a large component of the total dark matter. As reviewed in section [2], 
the axion can potentially also contribute to the overall dark matter content when it 
begins to oscillate at temperatures less than the reheating temperature of the saxion. 
This analysis only depends on the crude details of the saxion reheating temperature, 
which depending on the choice of inputs can sometimes be above, or below the axion 
temperature. As such, we will not dwell on this possibility in any great detail. 

In addition, although far less likely, it is also possible to consider scenarios where 
the saxion does not dominate the energy density. In such cases, the usual very strin- 
gent overclosure bounds on the gravitino relic abundance apply, effectively requiring 
Tp,fj < 10^ GeV for a gravitino of mass 7713/2 ~ 10 — 100 MeV. This can occur, for 
example, when the initial amplitude of the saxion is far smaller, at around the scale 
set by the decay constant of the axion. Nevertheless, because this requires a signif- 
icant calibration of various parameters such as the initial reheating temperature, in 
this paper we study the most natural situation where the saxion comes to dominate 
the energy density of the Universe. 

While the saxion neatly resolves some puzzling features typically found in grav- 
itino cosmology, it also has the potential to introduce additional issues. One such 
issue is that its decay products must not disrupt BBN. Indeed, as reviewed in section 
[2], the tight restrictions on the number of additional relativistic species which can 
be present translate into a constraint on the decay channels of the saxion. We find 
that this either requires introducing some additional weakly interacting particle into 
which the saxion can decay, or that the mass of the saxion must be sufficiently high 
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so that other decay channels to MSSM particles become kinematically available. 

Because the decay of the saxion indiscriminately dilutes various relics, it is im- 
portant to check whether an appropriate baryon asymmetry can be generated in F- 
theory GUT models. Rather than presenting an impediment, the decay of the saxion 
appears to open up more possibilities for generating a suitable baryon asymmetry! 
This is due to two key features. First, because the relic abundance of gravitinos can 
quite naturally fall in a viable range, the usually very tight prohibition on increas- 
ing the initial reheating temperature T^jj^ is no longer present. Rather importantly, 
many mechanisms for generating a large baryon asymmetry require a high value 
for T^H, which are commonly thought to be in conflict with the requirements of 
gravitino cosmology. Upon dispensing with the "gravitino problem", these options 
are once again available! Having said this, in typical models, the generated baryon 
asymmetry expected from such mechanisms is sometimes too large. In the present 
context, the same dilution effect discussed earlier turns out to retain an appropriate 
baryon asymmetry from scenarios such as standard leptogenesis in a natural range 
of parameters for F-theory GUT models. 

The rest of this section is organized as follows. Because it will play a very central 
role in much of the analysis to follow, we begin by analyzing the history of the saxion, 
and determine the precise conditions necessary for saxion domination. After this, we 
study the relic abundance of the gravitino. To this end, we first frame the discussion 
in general terms, asking under which circumstances we can expect the decay of a 
cosmological modulus to render the gravitino relic abundance independent of the 
temperatures T^^^, T^^^ and T^^. We find that there is a remarkably small range of 
parameters available, which are in fact consistent with the most naive expectations 
from F-theory GUTs! After this motivation, we present further details of how the 
saxion of F-theory GUTs satisfies all of these requirements. Next, we discuss the de- 
cay products generated by the decay of the saxion. After analyzing constraints from 
BBN, we show that standard leptogenesis can generate a suitable baryon asymmetry 
in the present class of models, with no fine tuning. 

4.1 Cosmology of the F-theory GUT Saxion 

Having outlined the general cosmology of F-theory GUTs, in this subsection we 
detail the primary role that the saxion plays as a cosmological modulus. In order for 
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the saxion to dilute the rehc abundance of a species such as the gravitino, it must 
be sufficiently long-lived in order for the coherent oscillations of the saxion field to 
dominate the energy density of the Universe. As obtained in section [3], the lifetime 
of the saxion is given by: 

T,a. - 2 X 10-6 sec ( ^"apJ"^ ) ' ^^^^^ 

Due to its long lifetime, the initial amplitude of the saxion field and its subsequent 
coherent oscillation can potentially dominate over other contributions to the energy 
density of the Universe. 

Because the saxion is a cosmological modulus, much of the general discussion 
reviewed in section [2] can now be applied to this special case of interest. Letting 
So denote the initial amplitude of the saxion, the temperature at which the saxion 
begins to oscillate is given by equation fl90l) so that: 



T:,, ~ 0.3 ■ ^/m^^^MpL = 0.6 ■ ^^pqMp^ ~ lO^^ GeV ■ (j^^^) 



1/2 

, (193) 



where in the second equality we have used the relation between the mass of the saxion 
and ApQ provided by equation (11701) . As anticipated, we note that this provides a 
link between cosmology and the PQ deformation. 

Once the saxion begins to oscillate, it will continue to do so until it decays. Using 
the general formula for the decay of a modulus, the associated temperature of decay 
is given by: 

/A \ 3/2 

T'rh = T^cay ~ 0.5 ■ v^rW^pI ~ 0.4 GeV ■ BjiH [j^^^j , (194) 

where by a small abuse of notation, we have identified the decay temperature with a 
"reheating temperature" , which is strictly speaking only correct if the saxion comes 
to dominate the energy density. Note in particular that T^^cay lower than T^^^. 
In addition, we also observe that the decay temperature of the saxion falls above 
the temperature required for BBN. Taking the maximal allowed branching ratio to 
axions consistent with BBN (which we will discuss later) so that Bg^aa ~ and 
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for a representative value of Apg ~ 100 GeV, we obtain the crude estimate!^ 

Tt,H ~ 1 GeV. (195) 

In between the temperature at which it begins to oscillate and decay, the saxion 
may come to dominate the energy density of the Universe. Letting T^^^ denote this 
temperature, this amounts to the condition: 

^osc ^ ^dom ^ ^decayi (196) 

where as reviewed in section [2] for a general cosmological modulus, the temperature 
Tdoni is given by: 

Tlom ~ S- min(T°^, T„^,J. (197) 

There are a priori two natural scales associated to the initial amplitude of sq. 
Because the saxion localizes on a matter curve with characteristic mass scale Mx ~ 
10^^-^ GeV, it is natural to take: 

So ~ Mx ~ 10^^-^ GeV. (198) 

Comparing our expressions for T^^^, TJ^„ and T^^cayy it follows that inequality fll96p 
indeed holds for sq ~ 10^^'^ GeV. We note in passing that from the perspective of 
the effective field theory, it is also possible to consider smaller field ranges set by the 
value of the axion decay constant so that sq ~ /a ~ 10^^ GeV. Note that in this 
case, TJ^m is smaller than TJ^^^^, indicating that the saxion in this case never comes 
to dominate the energy density of the Universe. 

Restricting to the most natural scenario where the saxion does indeed come to 
dominate the energy density, it will release a large amount of entropy when it decays. 
The associated dilution factor for any cosmological modulus again applies to the 



^Although the saxion is a cosmological modulus, it is interesting to note that in scenarios where 
it dominates the energy density of the Universe, it avoids the usual "moduli problem "because 
fa <C Mqut- Indeed, for a typical modulus of mass ~ rUsax, the analogous decay rate is of the 
form ~ m^/A^, where A is of the GUT, or Planck scale. As a consequence, the corresponding 
reheating temperature would then be much lower, jeopardizing BBN. 
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special case of the saxion, with the result: 



D - ^ . Jr^o V (199) 

Treating separately the two cases T^^^ > T^^ and T^^^ < T'^^, we therefore obtain 
the following expressions for the overall dilution factor: 



-1/2 1 


AOi° GeV\ 


f— >aa 1 
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-1/2 j 


ao^^-s GeV 


i— >aa 1 


^ So 



^T„,.<TO, - 10- ■ ^^-^) (t7^)- (201) 



100 GeV 



An important feature of these expressions is the overall dependence of the dilution 
factor on the initial reheating temperature, T^fj. Indeed, we note that when T^fj < 
T^sc^ the dilution factor becomes more potent as T^^j increases. This continues until 
T^fj ~ T^scy which point, the dilution factor ceases to depend on T]j^. This 
situation closely parallels the T^^j dependence of the gravitino, to which we shall 
shortly turn. 

Finally, it is also convenient to introduce the minimal reheating temperature 
required in order for saxion dilution to occur. This is given by the value of T^fj at 
which the dilution factor first equals one (the case of no dilution). Solving for T^^^ 
in the equality: 

-1- = J^min ~ 2~^rt)~ 



yields: 



„2 T-u 
■'O RH 



Ml 



Td^.„ = TL - "^T^H- (203) 



Using the explicit numerical values found in this section, it follows that the typical 
value of Tn is on the order of 10^ GeV. 



4.2 The Saxion-Gravitino Connection 

In the previous section we described the primary features of the saxion. The main 
point we found is that for typical values of F-theory GUT parameters, the oscillation 
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Figure 3: Schematic plot of the dilution factor D of the F-theory GUT saxion, 
the gravitino relic abundance in the absence of dilution, Vt^^^^h? and the net relic 
abundance after taking account of the dilution factor of the saxion as a function of 
the initial reheating temperature T^^. The graph of the dilution is given with respect 
to a scale distinct from that for the relic abundances. The plot depicts the special 
case where the freeze out temperature for the gravitino Ty^ ~ T*^^, the temperature 
at which the saxion begins to oscillate. When is below Tl)^,^, there is no dilution 
factor {D = 1). In the special case where T^^^ ~ T^^^, for all values of T^^j > Tf)^^^^ 
the total relic abundance of gravitinos is independent of T^^j. See figures H] and [5] for 
schematic plots of scenarios where 7^ ^osc- 
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Figure 4: Schematic plot of the dilution factor D of the F-theory GUT saxion, 
the gravitino relic abundance in the absence of dilution, ^l^y^h'^ and the net relic 
abundance after taking account of the dilution factor of the saxion as a function 
of the initial reheating temperature T^^- The graph of the dilution is given with 
respect to a scale distinct from that for the relic abundances. The plot depicts the 
case where the freeze out temperature for the gravitino T^^ > T^^^, the temperature 
at which the saxion begins to oscillate. When T^^j < T^^.^, there is no dilution 
factor. Note that in this case, the total relic abundance of gravitinos increases for 
values of the initial reheating temperature such that T^^^ < T^^j < tL^. 
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Figure 5: Schematic plot of the dilution factor D of the F-theory GUT saxion, 
the gravitino relic abundance in the absence of dilution, ^l^y^h'^ and the net relic 
abundance after taking account of the dilution factor of the saxion as a function 
of the initial reheating temperature T^^- The graph of the dilution is given with 
respect to a scale distinct from that for the relic abundances. The plot depicts the 
case where the freeze out temperature for the gravitino T^^ < T^^^, the temperature 
at which the saxion begins to oscillate. When T^^j < T^^.^, there is no dilution 
factor. Note that the total relic abundance of gravitinos decreases as T^^^ increases 
in the range T^^ < T^^ < T^^^. 
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of the saxion eventually comes to dominate the energy density of the Universe. The 
subsequent decay of the saxion will then dilute the relic abundance of all particle 
species. In this subsection we study the effects of this dilution on the relic abundance 
of gravitinos. We note that the idea of solving the gravitino problem due to a 
late decaying field has appeared for example in [42,43]. In this regard, the quite 
natural way that these ideas automatically appear in F-theory GUTs provides strong 
motivation for this type of cosmological scenario. 

Recall that in section [2] we reviewed the fact that for a gravitino of mass m^i2 ~ 
10 — 100 MeV, the resulting relic abundance would at first appear to violate the usual 
over closure bound: 

nl/^h^ < 0.1. (204) 

We find that the dilution of the saxion naturally resolves this issue because of a 
remarkable confluence of elements involving the masses of the gravitino and saxion, 
and the initial amplitude of the saxion. The basic point is that the relic abundance 
of gravitinos is governed by equation fl83|) : 

where here, we have included explicitly the dilution factor associated with the decay 
of the saxion. In addition, the gravitino freeze out temperature is given by equation 
(1731) as: 

T/ ~ 2 X 10- GeV ■ (^^) ' ( '-^] ' . (206) 

V 10 MeV/ \ nig J ^ ' 

Note that the overall relic abundance is given by the minimum of the freeze out 
temperature T^^^ and the initial reheating temperature T^^j. Indeed, recall that 
for m3/2 ~ 10 MeV and nig ~ 1 TeV, consistency with the overclosure bound would 
appear to require T^jj^ < 10^ GeV. Quite curiously, the amount of dilution associated 
with the saxion is also limited by the minimum of its oscillation temperature, and 
T^fj so that: 

D TT^^^rtTT- (207) 



si min(T„^,^, 7^° 



RH) 
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In other words, the thermally produced relic abundance of gravitinos is: 



^3/2^ ~ 2.7 X 10' 



si I min(T„^,„TO^) ) VlO^o GeV J [ m^/^ J \l TeV 



(208) 

While a priori, the oscillation temperature of the saxion and the gravitino mass are 
unrelated, in the context of F-theory GUTs, comparing the oscillation temperature 
Tq^^ of equation (11931) with the freeze out temperature of the gravitino provided by 
equation (12061) . we have: 

Tf ^ 2 X 10- GeV ■ (^^) ' ( '-^] ' . (210) 

VlOMeV/ \ rUg J ^ ^ 

In other words for typical values of 1713/2 ~ 10—100 MeV, nig ~ 1 TeV and ApQ ~ 100 
GeV, we obtain the remarkable relation: 

TL ~ tI,, (211) 

which should hold as an order of magnitude estimate. Returning to the gravitino 
relic abundance of equation (12081) . we therefore find: 

n„,h ~ 0.1 . (^^^j [j^] {-^) (iTw) ■ 



As a consequence, the overall relic abundance of gravitinos is independent of T/ 
T' and TO ! 



3/2' 

Using our expression for T^jj given by equation (I194p finally yields: 



07 / 10^^-^GeV y / ApQ y/V lO MeV \ f nig ^ 
n,/2h -0.07.5_/,^ J (^YoOGivJ [-^) \TTev) ' ^^^^^ 

which without any fine tuning satisfies the overclosure bound, and in certain circum- 
stances could account for the observed dark matter density! 

Turning the discussion around, the relation T^^^ ~ T!^2 niay be viewed as a type 
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of "resonance condition" which preferentially selects a window of values for the mass 
of the gravitino. Setting T^^^ ~ T^^^ and solving for 1713/2 in terms of Apg, we find: 

„,,.10MeV.(^)'"(^) (214) 

as an order of magnitude estimate. An intriguing consequence of this formula is that 
the mass of the gravitino is fairly insensitive to the value of the PQ deformation 
parameter. 

It is important to note that in the above estimate we have neglected potential 
temperature dependent corrections to the mass of the gravitino. In general, one 
would expect these corrections to decrease the mass of the gravitino above ~ 
lO*^'^ GeV. Our estimate for Ty^ ~ 10^° GeV is close to this scale, which suggests 
that there should be a mild decrease in the actual mass of the gravitino at this 
temperature. Thus, equation fl214p should be viewed as a lower estimate for the 
gravitino mass, which is reassuring in that 10 MeV is on the lower end of the mass 
range for the gravitino in F-theory GUT models. As a consequence, the coincidence 
window is roughly in the range 7713/2 ~ 10 — 100 MeV, which is intriguingly in the 
natural range expected for F-theory GUTs. 

4.2.1 F-theory and a Confluence of Parameters 

The result of the previous section suggests that some of the most distressing features 
of gravitino cosmology find a natural resolution in the context of F-theory GUTs. The 
crucial feature of this analysis is that although the gravitino relic abundance has non- 
trivial dependence on T^jj, this is exactly counterbalanced by the T^^j dependence 
of the dilution factor derived from the saxion. 

Given this remarkable conspiracy, it is natural to ask whether more general models 
with a late decaying modulus </> could potentially exhibit similar properties. In fact, 
as we now explain, the required relations between the mass of this modulus, its 
associated decay width, and the mass of the gravitino are only satisfied in a small 
window of values. 

The condition that the dilution factor exactly cancel the T^^j dependence of 
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^3/2^^ requires that the analogue of equation ( I21ip must hold: 

TL - 4/2- (215) 

Using the explicit expression for the oscillation temperature of a modulus as well as 
the value of T^j^ given by equation (I210p . we obtain: 

O.S^AV^^ ~ TL - T.I2 ~ 2 X 10^° GeV ■ (y^^)' , (216) 

where for simplicity, we have set rrig ~ 1 TeV. The mass of the gravitino is therefore 
given by: 

^3/2 / \i/4 



0-6- TT^T^ , (217) 



lOMeV V 100 GeV. 

In other words, a cosmological modulus with mass on the order of 100 — 1000 GeV 



will lead to a gravitino abundance which is independent of T'^^,T^^^ and T/,^ when 



the mass of the gravitino is in the range 10 — 100 MeV! 

We can also deduce further properties of the modulus by examining the overall 
relic abundance of gravitinos. Under conditions where equation fl215p is satisfied, 
it is enough to consider the special case where T'^^ > T^^^^T^j^. In this case, the 
constraint ^3/2/^^ < 0.1 translates to the condition 

{^) < 0.1, (218) 

where in the above, denotes the dilution factor of the cosmological modulus. 
Using the explicit expression for D^p reviewed in section El this becomes: 



fi'^ , Z,^ ^ ^ pji^ <0.1, (219) 

^'^ \ ms V2 keV/ " ' ^ ' 



where as before, denotes the decay rate of the modulus and 0o denotes the initial 
amplitude of this field. Parameterizing the decay rate as: 



1 rnl 



(220) 



647r A2 ' 

where A is some characteristic scale associated with the dynamics of the modulus, 
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we obtain the bound: 

which suggests a range of parameters similar to those of F-theory GUTs. 
4.3 Decay Products of the Saxion 

In the previous subsection we observed that the decay of the saxion and the associated 
release of entropy modifies the expected relic abundance of gravitinos, finding a 
remarkable confluence between the saxion oscillation temperature and the freeze out 
temperature for the gravitino. In addition we also found that the resulting relic 
abundance of thermally produced gravitinos is quite close to saturating the total 
amount of dark matter. 

But while the decay of the saxion effectively dilutes previously generated thermal 
relics, the end products of its decay can re-introduce another source for these same 
particles! Due to their overall longevity, decays to gravitinos and axions comprise the 
main source of additional potential relics. In this subsection we compute the relic 
abundance from such "non-thermally produced" gravitinos and axions. Whereas 
the amount of non-thermally produced axions is entirely negligible, depending on 
the mass of the gravitino, non-thermal production of gravitinos can also contribute 
towards the total dark matter^ 

For simplicity, in this subsection we restrict our attention to scenarios where the 
saxion comes to dominate the energy density of the Universe. Besides being the 
primary case of interest for F-theory GUTs, in the other more specialized case where 
the saxion is a subdominant component of the overall energy density, the resulting 
relic abundance will on general grounds generate a smaller component of the total 
matter content in comparison with a scenario with an era of saxion domination. 

In a scenario where the saxion dominates the energy density of the Universe, the 
fraction of the energy density transferred to the i^^ decay product is dictated by the 

''This is to be contrasted with for example, the result of [32], where in that case the production 
of gravitinos from the decay of the field responsible for supersymmctry breaking leads to essentially 
all of the gravitino relic abundance because the branching ratio to gravitinos is somewhat higher 
in the "sweet spot" scenario. 
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branching ratio: 

S,^,, = ^ (222) 

J- sax 

where in the above Ts^u the decay rate of the saxion to the z*'* species, and as before, 
Tsax denotes the total decay rate of the saxion. In this case, the overall yield of the 
i^^ species is given by: 

\rNT _ i,after _ -^fee/ore fT'sax,before _ -J ^ RH 

J- i — — -Ds~*ii ' — ^-Ds^ii ' ■ yZ,ZO) 

Safter ^after ^before ^ "^s 

where here, "before" and "after" are in reference to times close to the decay of the 
saxion. Note that as usual, Yi is constant as the Universe subsequently evolves. The 
non-thermally produced relic abundance from the z*^ species is therefore given by: 



VPc.O / \Pcfl J 2 nisax 

We now compute the value of VLf^h? in terms of the branching ratio Bg^aa- 
Using the explicit expression for Bg^ip^^^^^^^ in terms of Bg^aa given by equation 
fll84p . plugging in the present values of pc,o and sq reviewed in section [2] as well 
as the value of T^^^ obtained in equation (11941) . the non-thermally produced relic 
abundance of gravitinos and axions are respectively given by: 

»ZW ~ 0.9 . C-^) (225) 

^l^k^ - 5 X 10- . BIS^ (^) (^) . (226) 

In the natural range of masses for F-theory GUTs, it follows that the relic abundance 
from axions is completely negligible. On the other hand, the relic abundance of 
non-thermally produced gravitinos can potentially play a more prominent role. For 
example, using the representative values B ~ 10~^, 7713/2 ~ 20 MeV and ApQ ~ 100 
GeV, the resulting relic abundance of non-thermally produced gravitinos is ~ 0.05. 
On the other hand, when the mass of the gravitino is closer to the end range of 
F-theory GUT values at ~ 100 MeV, the relic abundance is at most ~ 10% of the 
dark matter, so that such gravitinos could comprise an additional component of the 
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dark matterl3 This is to be contrasted with the case of the "sweet spot" model of 
supersymmetry breaking [32,45], where the late decay of the field responsible for 
supersymmetry breaking generates most of the dark matter abundance. 

4.4 BBN and F-theory GUTs 

As reviewed in section [2], it is important to ascertain whether a given extension of 
the Standard Model disrupts the successful predictions of BBN. In this regard, we 
have already seen that the saxion can significantly alter the evolution of the Universe 
at temperatures above the start of BBN. In the most typical F-theory GUT scenario 
with an era of saxion domination, we find that the resulting reheating temperature is 
somewhat higher than Tbbn- Indeed, this imposes only the mild constraint ApQ > 1 
GeV. A far more significant constraint from BBN originates from the bound on 
the number of relativistic species. In the context of models with minimal matter 
content, this translates into the condition that the saxion must be allowed to decay 
to additional species beyond the axion, which in turn imposes a lower bound on the 
mass of the saxion. In terms of the PQ deformation, this amounts to the condition: 

ApQ > 60 GeV. (227) 

Note that this bound is indeed in accord with the crude expectation that the value 
of ApQ is most naturally near the weak scale. 

A late-decaying NLSP can also potentially alter the expected abundances of light 
elements generated by BBN. While a full study of BBN is beyond the scope of the 
present paper, a cursory inspection of the recent literature suggests that in compar- 
ison with standard cosmology, in a range of parameters for the gravitino and NLSP 
favored by F-theory GUTs, the expected abundance of light elements, such as the 
typically problematic '^Li is in somewhat better accord with observation. This is due 
to the fact that the late decay of the NLSP in such scenarios can alter the reaction 

^Due to the fact that non-thermaUy produced gravitinos are created at temperatures only a few 
orders of magnitude different from BBN, it is interesting to ask whether the resulting relics are 
indeed "cold" or "warm" at the time of matter recombination. In principle, "hot" dark matter can 
disrupt structure formation. As dark matter candidates, both gravitinos and axions are typically 
both sufficiently non-relativistic at the time of matter recombination to constitute cold dark matter 
candidates. We refer the interested reader to [44] for further discussion on the distinctions between 
hot, warm and cold dark matter. 



73 



rates of BBN. 



4.4.1 Decay Channels of the Saxion 

In subsection (12.7. ip we reviewed the general constraint on how many additional 
relativistic species can contribute such that the predictions of BBN remain in ac- 
cord with observation. In this subsection we consider the special case associated 
with the decay products of the saxion, closely following the quite general analysis 
for saxion decays presented in [31]. There are in principle two possibilities, depend- 
ing on whether or not the oscillation of the saxion comes to dominate the energy 
density of the Universe. As throughout, we shall specialize to the case where the 
saxion dominates the energy density of the Universe. Although the latter scenario 
is indeed more typical in the context of F-theory GUTs, for the sake completeness, 
in this subsection we study both possibilities in turn. Using the general constraint 
on the total number of allowed relativistic species reviewed in subsection 12.7.11 the 
contribution to the energy density from saxion generated axions at the time of BBN 
must satisfy the inequality: 

Pr / BBN 

Assuming the saxion comes to dominate the energy density of the Universe, the 
energy density stored in the saxion is directly transferred into the background radi- 
ation, as well as the axions. As such, we obtain the estimate: 

-) ~ Bs^aa, (229) 

Pr / BBN 

with Bs-,aa the branching ratio of the saxion to two axions. Inequality f l228p therefore 
provides a bound on the branching fraction to axions: 

Bs^aa < (230) 

In the minimal F-theory GUT, the primary decay channel at small Apq is given by 
the saxion to two axions. In order to satisfy this bound, the F-theory GUT must 
include some additional mode into which the saxion can decay. While it is certainly 
possible to posit the existence of some additional mode beyond the ones already 
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present in the MSSM, in a minimal scenario, the saxion will decay to MSSM degrees 
of freedom. As found in section [3l when available, the primary decay channel is 
for the saxion to decay to two Higgs bosons. This decay is kinematically allowed 
provided the mass of the saxion is sufficiently large: 



Using the expected mass of the Higgs of ~ 115 GeV, this translates into the following 
lower bound on the PQ deformation: 



We note that this crude kinematic constraint is fairly independent of the details of 
a particular F-theory GUT, such as the number of messengers, and in this regard 
is likely to be quite robust. Although it is also possible to derive a similar bound 
in the case where the saxion does not come to dominate the energy density of the 
Universe, in the context of F-theory GUTs, this is a far less likely scenario. We refer 
the interested reader to [31] for further details on this special case. We note that in 
general, the bound from the decay products of the saxion provides a somewhat tighter 
lower bound on Apg in comparison to the bound obtained from the requirement 
Trh > Tbbn- 

4.4.2 Comments on the Abundance of "^Li 

While the abundance of the light nuclei H~^, D'^, T~^, ^He~^~^, ^He~^^ expected 
from BBN are all in reasonable accord with observation, there is also some tension 
between the abundance of ^Lz expected based on the Standard Cosmology, and the 
observed abundance, which is typically a factor of 0.2 — 0.5 smaller. In fact, as 
mentioned in section [21 recent studies of the MSSM in scenarios with a gravitino 
in the mass range 10 — 100 MeV, with a bino or stau NLSP have recently been 
studied in [37] and even more recently in [35]. In this regard, it is interesting to 
note that the range of mass parameters expected in F-theory GUTs based on crude 
particle physics considerations are in rough accord with these studies. Although it 
is beyond the scope of this paper to address such detailed properties of BBN, it is 
quite encouraging that in [35], in the context of a gauge mediation scenario with a 



(231) 



ApQ > 60 GeV. 



(232) 
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gravitino of mass ~ 80 MeV and a bino NLSP of mass ~ 200 GeV that the resulting 
abundance of '^Li appears to be in better agreement with observation. It would be 
interesting to study this issue in greater detail. 



4.5 Baryon Asymmetry 

In the previous subsection we explained that the additional particle content of F- 
theory GUTs does not appear to disrupt the reaction rates necessary in BBN. At a 
more basic level, however, it is important to verify that the primary input of BBN, 
namely a sufficient baryon number asymmetry: 

^obs ^ UB-riT^ ^ ^ UB-riTs ^ ^ . ^ g X lO"!", (233) 

has in fact been generated! Here, Yg denotes the net yield of baryons. 

The creation of a sufficient baryon asymmetry is an especially acute problem 
in supersymmetric models with a gravitino LSP. Indeed, as reviewed in section [21 
in models without an era of moduli domination, it is quite common to lower the 
initial reheating temperature T'^jj to avoid overdosing the Universe from the thermal 
production of gravitinos. This in turn imposes strong constraints on the available 
mechanisms which can generate an appropriate baryon asymmetry. For example, in 
the context of models where the neutrinos of the Standard Model develop a small 
mass via the seesaw mechanism, heavy right-handed neutrinos couple to the Higgs 
up and lepton doublet through the superpotential term: 

W D Xi^KMNi, + Mi^^^N'j^Nl (234) 

where in the above, Mmaj denotes the Majorana mass of the right-handed neutrino, 
and is the Yukawa matrix in the neutrino sector, which for simplicity we shall take 
to be a 3 X 3 matrix. In the specific context of minimal F-theory GUTs which incor- 
porate a seesaw mechanism, simple estimates for the mass of the lightest Majorana 
mass give [2]: 

Ml ~ 3 X ioi2±i-5 GeV. (235) 

In standard leptogenesis, the subsequent decay of the right-handed neutrino to 
the Higgs and lepton doublet generates an overall lepton number density which is 
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converted via sphaleron processes to a baryon asymmetry. In order for this decay 
process to generate a sufficient baryon number asymmetry, however, the initial re- 
heating temperature must be greater than the Majorana mass: 

T'rh > Mmaj. (236) 

Indeed, for lower values of the initial reheating temperature, the decay products of 
the right-handed neutrinos are too dilute to generate the required baryon asymmetry. 

But while the decay of the saxion dilutes the relic abundance of thermally pro- 
duced gravitinos, it will indiscriminately also dilute any pre-existing baryon asym- 
metry! In this section we analyze whether a sufficient baryon asymmetry can be 
generated once the dilution effects of the saxion are taken into account, focussing on 
standard leptogenesis. We find that in the typical range of parameters for F-theory 
GUTs, standard leptogenesis typically generates a surplus baryon asymmetry which 
is diluted to acceptable values by the decay of the saxion. Similar studies on the 
compatibility of standard leptogenesis with a late decaying field which solves the 
gravitino problem have appeared, for example, in [42,43]. Again, we find it very 
reassuring that F-theory GUTs provide a natural setting for realizing such scenarios. 

Although we do not do so here, it is also possible to consider scenarios based 
on Dirac leptogenesis. In this case, an analogue of the seesaw mechanism generates 
small Dirac masses for the neutrinos, where the decay of the heavy particle associated 
with this "Dirac seesaw" generates light left- and right-handed neutrinos. Due to 
the difference in the efficiency of their interactions rates, this again can generate 
a lepton asymmetry, which is again converted to a baryon asymmetry. Insofar as 
standard leptogenesis can generate a viable level of baryon asymmetry, models with 
a similar range of parameters can also generate a sufficient baryon asymmetry in 
Dirac leptogenesis scenarios. While it would be interesting to study other alternative 
mechanisms for generating a large baryon asymmetry, it is beyond the scope of the 
present work to perform such an analysis. Indeed, the primary aim of this section is 
to demonstrate that in F-theory GUTs, standard mechanisms already generate an 
appropriate baryon asymmetry, without any additional assumptions. 



77 



4.5.1 Review of Standard Leptogenesis 



We now briefly outline the main points of standard leptogenesis [46], following the 
review [4 7]. In extensions of the Standard Model which generate a suitable Majorana 
mass term for the light neutrinos via the seesaw mechanism, the decay of heavy right- 
handed neutrinos into leptons and Higgses can generate a lepton asymmetry. This 
process satisfies the Sakharov conditions reviewed in section [2] because by construc- 
tion, the lepton number violation present in the Majorana mass term is converted 
to a baryon number violation via sphaleron processes. The violation of C and CP is 
somewhat more delicate, and at leading order originates from one loop contributions 
to the decay: 



in the obvious notation. The necessity of the one loop contribution for C and CP 
violation can be established for example, by appealing to the optical theorem. The 
amount of CP violation can be characterized in terms of the parameter ei, which 
roughly measures the overall complex phase of this one loop contribution. As re- 
viewed in [47], the net yield of leptons produced from this decay is: 



where in the context of supersymmetric models, g^,{MSSM) ~ 228.75, and k, is 
the "washout" factor which quantifies to what extent the decay of heavy neutrinos 
occurs out of equilibrium. The washout factor is given by integrating the Boltzmann 
equations, and in the range relevant to us is given by [47]: 



/ + h, 



(237) 




(238) 



0<r < 10:ft:(r) 



1 



(239) 



10 < r < 10^ : K(r 



2Vr2 + 9 
0.3 



(240) 



r(logr) 



0.8 



10^ <r: K(r) ~ (O.lr)^/^ ■ exp 




(241) 
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where in the above, r denotes a parameter associated with the efficiency of the 
reaction [47]: 

_-l/2 ^^ ^t^ . ,, 

(242) 



H{Mi) 1.7-327r Mi ' 

with Fat^ the decay rate for the hghtest of the heavy right-handed neutrinos. Sphaleron 
processes at high temperatures will automatically convert the lepton number asym- 
metry into a baryon asymmetry. The yield of net baryons Yb from this conversion 
is' 

The amount of CP violation in a given model depends on the size of the hierarchy 
in the right-handed neutrino masses. For example, in the case of an "extreme" 
hierarchy, M2/M1, M3/M1 > 10'^, the amount of CP violation is essentially dictated 
by a single complex number, and ei satisfies the Davidson-Iberra bound [48,49]: 

TT- 1-111^ 3a 5m -Ml /^..x 

hjxtreme Hierarchical : ei 5- ^ = ^di, 244 

' ' ~ 167r {Huf '''' ^ ^ 

where in the above, {Hu) denotes the vev of the Higgs up, and 6m = mmax — n^'mm ~ 
0.05 eV is the mass splitting in the light-neutrino sector. Further, a is an order one 
parameter associated with the value of the Yukawa couplings in the neutrino sector. 
Finally, in the context of the large tan/? = {Hy) / {Hd) scenarios studied in [3], the 
actual value of {H^) is given by: 

{Hu) = w sin /? ~ ~ 246 GeV. (245) 

Away from this extreme limit, the mass matrix for the right-handed neutrinos 
plays a more essential role, and there can be further sources of CP violation. For 
example, figure 1 of [50] illustrates that even when M2/Mi,M^/Mi ~ 10, the value 
of ei can deviate from e^i by one order of magnitude. More generally, in the case of 
less hierarchical masses, the resulting bound on ei is [49, 50] : 



Less Hierarchical : |ei| < max ( eoi, ) ~ ^{^)- (246) 



Ml 



Plugging in all numerical factors, the resulting upper bound on the baryon asym- 



79 



metry in these two situations is: 



Extreme Hierarchical : 77^ < 5 x 10 
Less Hierarchical : t]^"* < 10~^ ■ k, 




(248) 



where in the above, the superscript on the baryon asymmetry reflects the fact that 
the decay of a cosmological modulus could in principle dilute the total amount of 
baryon asymmetry generated. In general, we can expect an interpolation between 
the extreme, and less hierarchical scenarios. We now show that the expected range of 
Majorana masses in F-theory GUTs quite comfortably flts with the observed baryon 
asymmetry. 

4.5.2 Saxion Dilution and Standard Leptogenesis 

In the previous subsection we reviewed the main features of standard leptogenesis, 
focussing in particular on the distinction between hierarchical and non-hierarchical 
Majorana masses. In this subsection we analyze the effect of saxion dilution on 
the net baryon asymmetry. As reviewed near equation (12351) . in F-theory GUTs 
where the right-handed neutrinos have large Majorana masses, the natural mass 
scale associated to such flelds is roughly ~ 3 x lo'^^^^-^ GeV. In the present context 
then, standard leptogenesis is most natural in such cases when the initial reheating 
temperature T^^j is at or above this range of energy scales. 

We now determine whether the dilution factor from the saxion decay is sufficiently 
small to avoid overclosure from gravitinos, whilst at the same time, sufficiently large 
to avoid completely diluting the necessary baryon asymmetry generated by standard 
leptogenesis. In the natural range of parameters for F-theory GUTs, the Majorana 
mass of the right-handed neutrinos M^aj ^ 3 x IQ^^^'^^-^ GeV is greater than the 
freeze out temperature of the gravitino T^^^ ~ 10^'^ GeV. As a consequence, the relic 
abundance of thermally produced gravitinos is given by: 



In the range 7713/2 ~ 10 — 100 MeV, it follows that the dilution factor is bounded 




(249) 
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above by: 

, /10MeV\ , , 

D < 2 X 10~" ■ . 250 

V "^3/2 ) 

Multiplying the baryon asymmetry estimated in equations (12470 and (12481) by the 
dilution factor, the baryon asymmetry is therefore bounded above by: 



Extreme Hierarchical : ?7_b ^ 10 ■ k \ — — -^q — ) (251) 



^1012 Gev^ 

Less Hierarchical : 77^ < 2 x 10"^ ■ k, (252) 

where in the above estimate, we have used the fact that ra-^ji > 10 MeV. To proceed 
further, we now estimate the overall size of the wash out factor, k by determining 
the value of the parameter r in equation (I242p . Even without a complete theory of 
neutrino flavor, for our present purposes, it is enough to use the order of magnitude 
estimate for the Yukawa couplings: 

(A.At)n -4/^^-8x10-3 (253) 

Combined with the value of g^{MSSM) ~ 228.75, the resulting value of r is: 

6 X 10^2 QgV 



Ml 



2 X 10^^•^ (254) 



where in the final estimate we have plugged in the explicit value of Mi suggested by 
F-theory GUTs [2]. 

In the extreme hierarchical case, we see that a viable baryon asymmetry is only 
possible provided Mi > 10^^ GeV and k, an order one parameter. Returning to 
equation (12541) . it follows that in this range r is indeed quite small, so that equation 
(I239p implies k ~ 1/6. Plugging this value of k into (125 ip . the baryon asymmetry is 
bounded above by: 

Extreme Hierarchical : t]b ^ 2 x 10"^^ ■ ( f ^ ) . (255) 

V 10^2 GeV / 
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Generating the observed baryon asymmetry in this extreme case would then require: 

Extreme Hierarchical : Mi ~ 10^^ GeV, (256) 

which is remarkably close to the upper bound on Mi expected in F-theory GUTs. 

Next consider the more natural case for F-theory GUTs where the Majorana 
masses are not extremely hierarchical. In this case, inequality fl252p is saturated 
provided k, ~ 10-2 - 10'^. Returning to equations fl239|) - fl2iTD . this range of values 
requires k to be in the second range so that 10 < r < 10®. Evaluating k (r) at some 
representative values, we have: 

K(r~10)~10-2 (257) 
K(r ~ 100) ~ lO-l (258) 

In other words, in the range 10 < r < 100, leptogenesis with less hierarchical masses 
generates the requisite baryon asymmetry. In terms of the Majorana mass Mi, this 
corresponds to the range: 

Less Hierarchical : lO" GeV < Ah < 10^^ GeV, (259) 

which is in the expected range of Majorana masses estimated in [2]! 

To summarize, when the masses of the heavy neutrinos are not extremely hierar- 
chical, we find that within a natural window of values for Mi, the resulting baryon 
asymmetry matches with the observed value. This is due to the interplay between 
the dilution due to the saxion, and the natural range of Majorana masses expected 
in F-theory GUTs. Indeed, F-theory GUTs elegantly reconcile the apparent tension 
between standard leptogenesis and the "gravitino problem" . 

4.6 Messenger Relics 

In the above sections, we have seen that the cosmology of F-theory GUTs is re- 
markably insensitive to the initial reheating temperature of the Universe, T^h- ■'■^ 
the specific context of high scale gauge mediation scenarios, though, it is natural to 
ask whether the relic abundance of messengers will overdose the Universe. Indeed, if 
produced from the thermal bath, the large mass of these particles in F-theory GUTs, 
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Mmess ~ 10 GeV would lead to a very large relic abundance which even the decay 
of the saxion cannot dilute to an acceptable value. One possibility would be to take 
T^fj < Mmess ~ 10^^ GeV. In such a scenario the messengers would not have been 
produced by this initial condition assumption. However, this would not be attrac- 
tive in our scenario, because we have seen that essentially all the relevant physics is 
independent of T^^fj. Moreover this upper bound on the value of T^//, is potentially 
in conflict with leptogenesis. A more natural assumption, in line with the spirit of 
the present paper is to assume that the messenger fields can decay to some lighter 
fields, such as those present in the MSSM. We are currently investigating explicit 
models of F-theory GUTs which take this feature into account [19]. 



5 Future Directions 

In this paper we have found that in F-theory GUTs, the gravitino and in some cases 
the axion can provide a prominent component of the total dark matter, which quite 
remarkably is independent of the initial reheating temperature of the Universe, T^^j. 
On the other hand, such candidates leave open the issue of accounting for the recent 
experimental results such as PAMELA [51], which could potentially be explained in 
terms of dark matter physics. In this regard, it is important to investigate whether 
F-theory GUTs provide additional dark matter candidates [19]. 

In the context of F-theory GUTs, the decay of the saxion which dilutes the 
abundance of gravitinos will also dilute the relic abundance of any other dark matter 
candidate. Assuming that the dark matter originates from a cold thermal relic, its 
abundance scales inversely with its cross section: 

^^DM/i' OC i -. (260) 

Letting M denote a characteristic mass scale associated with the dark matter, ctdm ~ 
M~^, which would naturally suggest a mass scale of order M ~ 1 TeV. Taking into 
account the dilution factor D ~ 10~^ — 10~^ from the decay of the saxion would 
instead suggest that M ~ 100 TeV. The downside to this is that this also lowers the 
cross section by a factor of 10~^ — 10~^. Unless there is a substantial enhancement 
either in its density or cross section at small velocities, this type of dark matter 
candidate would then be too small to be detectable in current dark matter searches. 
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Nevertheless, an exact analysis of potential dark matter candidates depends on 
the details of a given model. In this regard, the tight structure of F-theory GUTs 
also provides additional candidates associated with degrees of freedom located near 
the F-theory GUT seven-brane. For example, four- dimensional GUTs always include 
a U{1)b-l gauge boson. We are currently investigating the details of such a scenario 
[19]. 
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